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Abstract 

We improve the theoretical estimates of the critical exponents for the three- 
dimensional XY universality class. We find a = —0.0146(8), 7 = 1.3177(5), 
ly = 0.67155(27), t/ = 0.0380(4), (3 = 0.3485(2), and 5 = 4.780(2). We observe 
a discrepancy with the most recent experimental estimate of a; this discrep- 
ancy calls for further theoretical and experimental investigations. Our results 
are obtained by combining Monte Carlo simulations based on finite-size scal- 
ing methods, and high-temperature expansions. Two improved models (with 
suppressed leading scaling corrections) are selected by Monte Carlo computa- 
tion. The critical exponents are computed from high-temperature expansions 
specialized to these improved models. By the same technique we determine 
the coefficients of the small-magnetization expansion of the equation of state. 
This expansion is extended analytically by means of approximate parametric 
representations, obtaining the equation of state in the whole critical region. 
We also determine the specific-heat amplitude ratio. 

PACS Numbers: 05.70.Jk, 64.60.Pr, 75.10.Hk, 11.15.Me 
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I. INTRODUCTION 



In the theory of critical phenomena continuous phase transitions can be classified into 
universality classes determined only by a few properties characterizing the system, such as 
the space dimensionality, the range of interaction, the number of components of the order 
parameter, and the symmetry. Renormalization-group (RG) theory predicts that, within a 
given universality class, critical exponents and scaling functions are identical for all systems. 
Here we consider the three-dimensional XY universality class, which is characterized by a 
two-component order parameter, 0(2) symmetry, and short-range interactions. 

The superfluid transition of ^He, whose order parameter is related to the complex quan- 
tum amplitude of the helium atoms, belongs to the three-dimensional XY universality class. 
It provides an exceptional opportunity for an experimental test of the RG predictions, es- 
sentially because of the weakness of the singularity in the compressibility of the fiuid, of 
the purity of the samples, and of the possibility of performing the experiments, such as the 
Space Shuttle experiment reported in in a microgravity environment, thereby reducing 
the gravity- induced broadening of the transition. Because of these favorable conditions, 
the specific heat of liquid helium was accurately measured to within a few nK from the A 
transition, i.e., very deep in the critical region, where the scaling corrections to the expected 
power-law behavior are small. The experimental low-temperature data for the specific heat 
were analyzed assuming the behavior for t = (T — T^/T^ — * 



This result represents a challenge for theorists because its uncertainty is substantially smaller 
than those of the theoretical calculations. We mention the best available theoretical esti- 
mates of a: a = —0.0150(17) obtained using high-temperature (HT) expansion techniques 
0], a = —0.0169(33) from Monte Carlo (MC) simulations using finite-size scaling (FSS) 
techniques 0], and a = —0.011(4) from field theory [g]. 

The aim of this paper is to substantially improve the precision of the theoretical estimates 
of the critical exponents, reaching an accuracy comparable with the experimental one. For 
this purpose, we will consider what we call "improved" models. They are characterized by 
the fact that the leading correction to scaling is absent in the expansion of any observable 



^This value of A is close to the best available theoretical estimates, i.e., A = 0.529(8) from field 
theory Q and A = 0.531(14) from Monte Carlo simulations 

2 Ref. Q reported a = -0.01285(38) and A+/A~ = 1.054(1). But, as mentioned in footnote [15] of 
Ref. Q, the original analysis was slightly in error. Ref. Q reports the new estimates a = —0.01056 
and A+/j4^ = 1.0442. J. A. Lipa kindly communicated us ||5| that the error on a is the same as in 



Cnit) = (1 + C\t\^ + Dt)+B 

with A = l/2.[] This provided the estimate 



(1) 




(2) 



Ref. i. 
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near the critical point. Moreover, we will combine MC simulations and analyses of HT series. 
We exploit the effectiveness of MC simulations to determine, by using FSS techniques, the 
critical temperature and the parameters of the improved Hamiltonians, and the effectiveness 
of HT methods to determine the critical exponents for improved models, especially when 
a precise estimate of jSc is available. Such a combination of lattice techniques allows us to 
substantially improve earlier theoretical estimates. We indeed obtain 

a = -0.0146(8), (3) 

where, as we will show, the error estimate should be rather conservative. The theoretical 
uncertainty has been substantially reduced. We observe a disagreement with the experi- 
mental value (^. The point to be clarified is whether this disagreement is significant, or it 
is due to an underestimate of the errors reported by us and/or in the experimental papers. 
We think that this discrepancy calls for further theoretical and experimental investigations. 
A new-generation experiment in microgravity environment is currently in preparation 0; it 
should clarify the issue from the experimental side. 

In numerical (HT or MC) determinations of critical quantities, nonanalytic corrections 
to the leading scaling behavior represent one of the major sources of systematic errors. 
Considering, for instance, the magnetic susceptibility, we have 

X = Cr'' (1 + ao,it + ao,2t^ + ••• + + 01,2^^^ + ... -t- 02,1^^' + ...) . (4) 

The leading exponent 7 and the correction-to-scaling exponents A,A2,..., are universal, 
while the amplitudes C and nonuniversal. For three-dimensional XY systems, the 

value of the leading correction-to-scaling exponent is A ^ 0.53 0,0, and the value of the 
subleading exponent is A2 ~ 2 A 0. 

The leading nonanalytic correction is the dominant source of systematic errors in 
MC and HT studies. Indeed, in MC simulations the presence of this slowly-decreasing term 
requires careful extrapolations, increasing the errors in the final estimates. In HT studies, 
nonanalytic corrections introduce large and dangerously undetectable systematic deviations 



in the results of the analyses. Integral approximants (see, e.g., Ref. []TOl for a review) 
can in principle cope with an asymptotic behavior of the form (^; however, in practice, 
they are not very effective when applied to the series of moderate length available today. 
Analyses meant to effectively allow for the leading confluent corrections are based on biased 
approximants, where the value of j3c and the first non-analytic exponent A are introduced as 
external inputs (see e.g. Refs. [|TT]-|T^). Nonetheless, their precision is still not comparable 



to that of the experimental result (0), see e.g. Ref. |jT3|. The use of improved Hamiltonians, 
i.e., models for which the leading correction to scaling vanishes (ai i = in Eq. (§)|), can 
lead to an additional improvement of the precision, even without a substantial extension of 
the HT series. 

The use of improved Hamiltonians was first suggested in the early 80s by Chen, Fisher, 
and Nickel jl^ who determined improved Hamiltonians in the Ising universality class. The 



^Actually, for improved models, ai_j = for all is. 
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crux of the method is a precise determination of the optimal value of the parameter appearing 
in the Hamiltonian. One can determine it from the analysis of HT series, but in this case it 
is obtained with a relatively large error 



and the final results do not significantly 
improve the estimates obtained from standard analyses using biased approximants. 

Recently pT|-p^,p|,p^ it has been realized that FSS MC simulations are very effective 
in determining the optimal value of the parameter, obtaining precise estimates for several 
models in the Ising and XY universality classes. The same holds true of models in the 0(3) 
and 0(4) universality classes p5[. Correspondingly, the analysis of FSS results obtained 
in these simulations has provided significantly more precise estimates of critical exponents. 
An additional improvement of the precision of the results has been obtained by combining 
improved Hamiltonians and HT methods. Indeed, we already showed that the analysis of HT 
series for improved models [p!^,|6|,p6| provides estimates that are substantially more precise 
than those obtained from the extrapolation of the MC data alone. 

In this paper we consider again the XY case. The progress with respect to the studies 
of Refs. [|,^] is essentially due to the improved knowledge of /3c and of the parameters of 
the improved Hamiltonians obtained by means of a large-scale MC simulation. The use of 
this information in the analysis of the improved HT (IHT) series allows us to substantially 
increase the precision and the reliability of the results, especially of the critical exponents. 
As we shall see, in order to determine the critical exponents, the extrapolation to (3c of 
the IHT series, using biased integral approximants, is more effective than the extrapolation 
L oo of the FSS MC data. Moreover, we consider two improved Hamiltonians. The 
comparison of the results from these two models provides a check of the errors we quote. 
The estimates obtained for the two models are in very good agreement, providing support for 
our error estimates and thus confirming our claim that the systematic error due to confluent 
singularities is largely reduced when analyzing IHT expansions. 

We consider a simple cubic (sc) lattice, two-component vector fields (j)x = {(l>x \ and 
two classes of models depending on an irrelevant parameter: the 0^ lattice model and the 
dynamically diluted XY (dd-XY) model. 

The Hamiltonian of the 0^ lattice model is given by 



The dd-XY model is defined by the Hamiltonian 



n 



•dd 



(5) 



(6) 



by the local measure 



and the partition function 



27r 



(7) 



(8) 
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In the limit D oo the standard XY lattice model is recovered. We expect the phase tran- 
sition to become of first order for D < Dtri- -Dtri vanishes in the mean-field approximation, 
while an improved mean-field calculation based on the "star approximation" of Ref. ||27 
gives Dtj-i < 0, so that we expect D^ri < 0. 

The parameters A in 7i^4 and D in Tidd can be tuned to obtain improved Hamiltonians. 
We performed an accurate numerical study, which provided estimates of A*, D*, of the inverse 
critical temperature jSc for several values of A and D, as well as estimates of the critical 
exponents. Using the linked-cluster expansion technique, we computed HT expansions of 
several quantities for the two theories. We analyzed them using the MC results for A*, D* 
and Pc, obtaining very accurate results, e.g., Eq. (^). 

We mention that the 0^ lattice model Ti^i has already been considered in MC and HT 
studies P,P,p6| . With respect to those works, we have performed additional MC simulations 
to improve the estimate of A* and determine the values of f3c- Moreover, we present a new 
analysis of the IHT series that uses the MC estimates of Pc to bias the approximants, leading 
to a substantial improvement of the results. 

In Table | we report our results for the critical exponents, i.e., our best estimates obtained 
by combining MC and IHT techniques — they are denoted by MC+IHT — together with the 
results obtained from the analysis of the MC data alone. There, we also compare them 
with the most precise experimental and theoretical estimates that have been obtained in the 
latest years. When only z/ or a is reported, we used the hyperscaling relation 2 — a = Su to 
obtain the missing exponent. Analogously, if only 77 or 7 is quoted, the second exponent was 
obtained using the scaling relation 7 = (2 — 77)//; in this case the uncertainty was obtained 
using the independent-error formula. The results we quote have been obtained from the 
analysis of the HT series of the XY model (HT), by Monte Carlo simulations (MC) or 
by field-theory methods (FT). The HT results of Ref. have been obtained analyzing 
the 21st-order HT expansions for the standard XY model on the sc and the bcc lattice, 
using biased approximants and taking j3c and A from other approaches, such as MC and 
FT. The FT results of Refs. P,p8| have been derived by resumming the known terms of 



the fixed-dimension g expansion: the /3 function is known to six-loop order |29|, while the 
critical-exponent series are known to seven loops ||30|. The estimates from the e expansion 
have been obtained resumming the available O(e^) series |l3T| , |32[| . 



We also present a detailed study of the equation of state. We first consider its expan- 
sion in terms of the magnetization in the high-temperature phase. The coefficients of this 
expansion are directly related to the zero-momentum n-point renormalized couplings, which 
were determined by analyzing their IHT expansion. These results are used to construct 
parametric representations of the critical equation of state which are valid in the whole 
critical region, satisfy the correct analytic properties (Griffiths' analyticity) , and take into 
account the Goldstone singularities at the coexistence curve. From our approximate repre- 
sentations of the equation of state we derive estimates of several universal amplitude ratios. 
The specific-heat amplitude ratio is particularly interesting since it can be compared with 
experimental results. We obtain A'^/A^ = 1.062(4), which is not in agreement with the 
experimental result A'^ /A~ = 1.0442 of Refs. [00. It is easy to trace the origin of the dis- 
crepancy. In our method as well as in the analysis of the experimental data, the estimate of 
A'^/A~ is strongly correlated with the estimate of a. Therefore, the discrepancy we observe 
for this ratio is a direct consequence of the difference in the estimates of a. 



5 



TABLE I. Estimates of the critical exponents. See text for the explanation of the symbols in 
the second column. We indicate with an asterisk (*) the estimates that have been obtained using 
the hyperscaling relation 2 — a = or the scaling relation 7 = (2 — r])^. 



Ref. 


Method 


7 


V 




a 


this work 


MC+IHT 


1.3177(5) 


0.67155(27) 


0.0380(4) 


-0.0146(8)* 


this work 


MC 


1.3177(10)* 


0.6716(5) 


0.0380(5) 


-0.0148(15)* 


i 


(2000) 


IHT 


1.3179(11) 


0.67166(55) 


0.0381(3) 


-0.0150(17)* 


P 


(1999) 


HT 








-0.014(9), -0.022(6) 


11 


(1997) 


HT, sc 


1.325(3) 


0.675(2) 


0.037(7)* 


-0.025(6)* 






HT, bcc 


1.322(3) 


0.674(2) 


0.039(7)* 


-0.022(6)* 


i 


(1999) 


MC 


1.3190(24)* 


0.6723(11) 


0.0381(4) 


-0.0169(33)* 


M 


(1999) 


MC 


1.315(12)* 


0.6693(58) 


0.035(5) 


-0.008(17)* 


11 

LL- II 


(1996 


MC 


1.316(3)* 


0.6721(13) 


0.0424(25) 


-0.0163(39)* 


ii 


(1995) 


MC 




0.6724(17) 




-0.017(5)* 




(1993) 


MC 


1.307(14)* 


0.662(7) 


0.026(6) 


-0.014(21)* 




(1990) 


MC 


1.316(5) 


0.670(2) 


0.036(14)* 


-0.010(6)* 


M 


(1999) 


FT d = 3 exp 


1.3164(8) 


0.6704(7) 


0.0349(8) 


-0.0112(21) 


i 


(1998) 


FT d = 3 exp 


1.3169(20) 


0.6703(15) 


0.0354(25) 


-0.011(4) 


i 


(1998) 


FT e-exp 


1.3110(70) 


0.6680(35) 


0.0380(50) 


-0.004(11) 


ii 


(1996) 


^He 




0.67019(13)* 




-0.01056(38) 


m 


(1993) 


^He 




0.6705(6) 




-0.0115(18)* 


m 


(1992) 


^He 




0.6708(4) 




-0.0124(12)* 


m 


(1984) 


^He 




0.6717(4) 




-0.0151(12)* 


m 


(1983) 


^He 




0.6709(9)* 




-0.0127(26) 



Finally, we also discuss the two-point function of the order parameter, i.e., the structure 
factor, which is relevant in scattering experiments with magnetic materials. 

The paper is organized as follows. In Sec. we present our Monte Carlo results. After 
reviewing the basic RG ideas behind our methods, we present a determination of the im- 
proved Hamiltonians and of the critical exponents. We discuss the several possible sources 
of systematic errors, and show that the approximate improved models we use have signifi- 
cantly smaller corrections than the standard XY model. A careful analysis shows that the 
leading scaling corrections are reduced at least by a factor of 20. We also compute to high 
precision for several values of A and D\ this is an important ingredient in our IHT analyses. 
Details on the algorithm appear in App. 0. 

In Sec. |TTl| we present our results for the critical exponents obtained from the analysis 
of the IHT series. The equation of state is discussed in Sec. After reviewing the basic 
definitions and properties, we present the coefficients of the small-magnetization expan- 
sion, again computed from IHT series. We discuss parametric representations that provide 
approximations of the equation of state in the whole critical region and compute several 
universal amplitude ratios. In Sec. |V| we analyze the two-point function of the order param- 
eter. Details of the IHT analyses are reported in App. p. The definitions of the amplitude 
ratios we compute can be found in App. |C[ 
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II. MONTE CARLO SIMULATIONS 



A. The lattice and the quantities that were measured 



We simulated sc lattices of size V = L^, with periodic boundary conditions in all three 
directions. In addition to elementary quantities like the energy, the magnetization, the spe- 
cific heat or the magnetic susceptibility we computed so-called phenomenological couplings, 
i.e., quantities that, in the critical limit, are invariant under RG transformations. They are 
well suited to locate the inverse critical temperature Pc- They also play a crucial role in the 
determination of the improved Hamiltonians. 

In the present study we consider four phenomenological couplings. We use the Binder 
cumulant 



^4: — / ^2\2 ' 

and the analogous quantity with the 6th power of the magnetization 

((m2)3) 



2\3 



(9) 



(10) 



where m = ^ (px is the magnetization of the system. 

We also consider the second-moment correlation length divided by the linear extension 
of the lattice ^2nd/L. The second-moment correlation length is defined by 



^2nd 



X/F~l 
4sin(7r/L)2' 



where 



X = 




(11) 



(12) 



is the magnetic susceptibility and 

1 



F = 



V 



exp 



271X1 



(13) 



is the Fourier transform of the correlation function at the lowest non-zero momentum. 

The list is completed by the ratio Za/Zp of the partition function Za of a system with 
anti-periodic boundary conditions in one of the three directions and the partition function 
Zp of a system with periodic boundary conditions in all directions. Anti-periodic boundary 
conditions in the first direction are obtained by changing sign to the term 4'x4'y of the 
Hamiltonian for links (xy) that connect the boundaries, i.e., for x = {L,X2,X3) and y = 
(l,a;2,a;3). The ratio Za/Zp can be measured by using the boundary-flip algorithm, which 
was applied to the three-dimensional Ising model in Ref. and generalized to the XY 
model in Ref. |Q. As in Ref. in the present work we used a version of the algorithm 



that avoids the flip to anti-periodic boundary conditions. For a detailed discussion see App. 



7 



B. Summary of finite-size methods 



In this subsection we discuss the FSS methods we used to compute the inverse critical 
temperature, the couphngs A* and D* at which leading corrections to scaling vanish, and 
the critical exponents v and rj. 



1. Summary of basic RG results 

The following discussion of FSS is based on the RG theory of critical phenomena. We 
first summarize some basic results. In the three-dimensional XY universality class there exist 
two relevant scaling fields Ut and m/j, associated to the temperature and the applied field 
respectively, with RG exponents yt and yh- Moreover, there are several irrelevant scaling 
fields that we denote by -Uj, i > 3, with RG exponents > 7/3 > ?/4 > 1/5 > ... . 

The RG exponent 7/3 = —uj of the leading irrelevant scaling field M3 has been com- 
puted by various methods. The analysis of field-theoretical perturbative expansions 
gives oj = 0.802(18) (e expansion) and u = 0.789(11) {d = 3 expansion). In the present 
work we find a result for uj that is consistent with, although less accurate than, the field- 
theoretical predictions. We also mention the estimate u = 0.85(7) that was obtained by 
the "scaling-field" method, a particular implementation of Wilson's "exact" renormalization 
group. Although it provides an estimate for u that is less precise than those obtained from 
perturbative field-theoretic methods, it has the advantage of giving predictions for the ir- 
relevant RG exponents beyond y^. Ref. predicts 1/4 = —1.77(7) and y^ = —1.79(7) (1/421 
and 7/422 in their notation) for the XY universality class. Note that, at present, there is no 
independent check of these results. Certainly it would be worthwhile to perform a Monte 
Carlo renormalization group study. With the computational power available today, it might 
be feasible to resolve subleading correction exponents with a high-statistics simulation. 

In the case of U4, Uq, and ^2nd/ L we expect a correction caused by the analytic background 
of the magnetic susceptibility. This should lead to corrections with y^ = —{2 — 7]) ^ —1.962. 
We also expect corrections due to the violation of rotational invariance by the lattice. For 
the XY universality class, Ref. |^ predicts y-j = —2.02(1). Note that the numerical values 



of ye and 7/7 are virtually identical and should hence be indistinguishable in the analysis of 
our numerical data. 

We wish now to discuss the FSS behavior of a phenomenological coupling R; in the 
standard RG framework, we can write it as a function of the thermal scaling field Ut and of 
the irrelevant scaling fields Ui. For L — > 00 and (3 — >■ /3c(A), we have 

R{L, (3, A) = roiutLy^) + ^ r,(MiL^') + ... , (14) 

j>3 

where we have neglected terms that are quadratic in the scaling fields of the irrelevant 
operators, i.e., corrections of order L^^^ ^ L~^'^. Note that we include here the corrections 
due to the analytic background (with exponent L~y^ ^ L''"^). In the case of f/4, f/g, and 
^2nd/L (but not Za/Zp), in Eq. (|14D we have also discarded terms of order ^ L~^-^ . 

The functions tq^z) and ri(z) are smooth and finite for 2; — 0, while Mt(/3, A) and Ui{[3, A) 
are smooth functions of [3 and A. Note that, by definition, Mf(/3, A) ~ 13 — l3c{X)- In the limit 
t ^ and utiy^ ~ (/3 - (30)1^'" 0, we can further expand Eq. (P^, obtaining 
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R{L, (3, X) = R* + ct{(3, A) L^' + J] c,{(3, A) L^' + 0{{(3 - (3,)^L^y\ L^y^ ,11^^+^'), (15) 

i 

where R* = ro(0) is the value at the critical point of the phenomenological coupling. 

2. Locating f3c 

We locate the inverse critical temperature I3c by using Binder's cumulant crossing method. 
This method can be applied in conjunction with any of the four phenomenological couplings 
that we computed. 

In its simplest version, one considers a phenomenological coupling R{(3, L) for two lattice 
sizes L and L' = bL. The intersection /Across of the two curves R{P, L) and R{l3, L') provides 
an estimate of jSc- The convergence rate of this estimate /Across towards the true value can 
be computed in the RG framework. 

By definition, /Across at fixed b, L, and A is given by the solution of the equation 

R{L,(3,X) = R{bL,(3,X). (16) 

Using Eq. (|I4D, one immediately verifies that /Across converges to Pc faster than L^^*. Thus, 
for L — i> oo, we can use Eq. (ITB) and rewrite Eq. (ITB]) as 



CtiP, A) L^^ + ciP, A) ^ CtiP, A) {bL)y^ + ciP, A) (6L)^^ (17) 

Then, we approximate q(/5. A) ~ (/? — Pc) and Cj(/5, A) ~ Cj(/?c, A) = Cj. Remember that 
Ct{Pc, A) = by definition. Using these approximations we can explicitly solve Eq. (|17[) with 
respect to P, obtaining 

C3 (1 - by^) Ly^ , , 

The leading corrections vanish like L~y^^y'-^ ^ L^^'^. Inserting /Across into Eq. (0), we obtain 

i?cross = R*+ ly^ + ... , (19) 

which shows that the leading corrections vanish like L^^. 

Given a precise estimate of R*, one can locate Pc from simulations of a single lattice size, 
solving 

R{L,p) = R\ (20) 
where the corrections vanish like L~y^^y^. 
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3. Locating A* and D* 

In order to compute the value A* for which the leading corrections to scaling vanish, we 
use two phenomenological couplings Ri and R2. First, we define l3f{L, A) by 

i?i(L,/3;,A) = i?ij, (21) 

where Rij is a fixed value, which we can choose freely within the appropriate range. It is 



easy to see that /5/(L, A) — *• /9c(A) as L — 00. Indeed, using Eq. (p!4D, we have 



/3^(L, A) - /3,(A) = z^L~y^ - !mM«3(/3,)L^3-.. + ... ^ (22) 

where we have used yt > lysl and Zf is defined as the solution of ri^o{zf) = Rij- We have 
added a subscript 1 to make explicit that all scaling functions refer to Ri. If Rij ^ R*, we 
can expand the previous formula, obtaining 

Pf{L, A) ^ /3,(A) + ^hl^L-y^ _ ^^-y*+.3. (23) 

Notice that for Rij = RI the convergence is faster, and thus we will always take Ri j ^ RI. 
Next we define 

R{L,X)=R2{L,(]f,X). (24) 
For L ^ 00 and Rij ^ RI, we have 

Rm X)=R; + ^(i?l,/ -Rl)+J2 f C2,i - ^ C; 

= i?* + 5^Q(A)L^% (25) 

i 

which shows that the rate of convergence is determined by L^^. 

In order to find A*, we need to compute the value of A for which Cj(A) = 0. We can 
obtain approximate estimates of A* by solving the equation 

R{L,X) = R{bL,X). (26) 

Using the approximation (^) one finds 

Across = A*-^^4iL" + ..., (27) 

where Cg is the derivative of C3 with respect to A, and 

Ky _ 

R,ro.. = R* -CA^-^Ly\ (28) 
I -hi 

In principle, any pair R2 of phenomenological couplings can be used in this analysis. 
However, in practice we wish to see a good signal for the corrections. This means, in 
particular, that in c^^ 02,3 — c^ ^ Ci_3 the two terms should add up rather than cancel. Of 
course, also the corrections due to the subleading scaling fields should be small. 
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4- The critical exponents 



Typically, the thermal RG exponent j/f = 1/z/ is computed from the FSS of the derivative 
of a phenomenological couphng R with respect to /? at Pc- Using Eq. (|1^) one obtains 



dR 
'dp 



=3 i=3 



Hence, the leading corrections scale with L^^. However, in improved models in which 
^3(/3c) = 0, the leading correction is of order L^*. Note that corrections proportional to 
l^y-i-yt ~ /,-2.3 g^jjj present even if the model is improved. In Ref. |2l|, for the spin-1 



Ising model, an effort was made to eliminate also this correction by taking the derivative 
with respect to an optimal linear combination of f3 and D instead of (3. Here we make no 
attempt in this direction, since corrections of order L~^'^ are subleading with respect to 
those of order L^* ^ 

In practice it is difficult to compute the derivative at (3c, since f3c is only known numer- 
ically, and therefore, it is more convenient to evaluate OR/ 3(3 at /9/ (see Eq. (|2l|) ). This 
procedure has been used before, e.g., in Ref. [^. In this case, Eq. ( pi3| ) still holds, although 



with different amplitudes that depend on the particular choice of the value of i?/. 

The exponent r] is computed from the finite-size behavior of the magnetic susceptibility, 

i.e., 

X\p, oc L'-\ (30) 

Also here the corrections are of order L^^ for generic models, and of order L^'* for improved 
ones. 



5. Estimating errors caused by residual leading scaling corrections 



In Ref. [^, the authors pointed out that with the method discussed in Sec. [11 B 3| , the 
leading corrections are only approximately eliminated, so that there is still a small leading 
scaling correction which causes a systematic error in the estimates of, e.g., the critical 
exponents. The most naive solution to this problem consists in adding a term L~'^ to the 
fit ansatz, i.e., in considering 



dR 

'8(3 



(31) 



However, by adding such a correction term, the precision of the result decreases, so that 
there is little advantage in using (approximately) improved models. A more sophisticated 
approach is based on the fact that ratios of leading correction amplitudes are universal |^ . 

Let us consider a second phenomenological coupling R = R{L, (3f, A), which, for L oo, 
behaves as 



R = R* + C3L- 



(32) 
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The universality of the correction amphtudes imphes that the ratio B/c^ is the same for 
the 0*^ model and the dd-XY model and is independent of A and D. Therefore, this ratio 
can be computed in models that have large corrections to scaling, e.g., in the standard XY 
model. Then, we can compute a bound on B for the (approximately) improved model from 
the known ratio B/c^ and a bound for C3. This procedure was proposed in Ref. [pT| . 



C. The simulations 

We simulated the 0"^ and the dd-XY model using the wall-cluster update algorithm of 
Ref. combined with a local update scheme. The update algorithm is discussed in detail 
in App. 1^, where we also report an analysis of its performance and the checks we have done. 

Most of the analyses need the quantities as functions of (3. Given the large statistics, we 
could not store all individual measurements of the observables. Therefore, we did not use 
the reweighting method. Instead, we determined the Taylor coefficients of all quantities of 
interest up to the third order in {(3 — Ps), where Ps is the value of /? at which the simulation 
was performed. We checked carefully that this is sufficient for our purpose. For details, see 
App. 0. 

Most of our simulations were performed at A = 2.1 in the case of the 0^ model and at 
D = 1.03 in the case of the dd-XY model. A = 2.1 is the estimate of A* of Ref. 0, and 
D = 1.03 is the result for D* of a preliminary analysis of MC data obtained on small lattices. 

In addition, we performed simulations at A = 2.0 and 2.2 for the 0^ model and D = 0.9 
and 1.2 for the dd-XY model in order to obtain an estimate of the derivative of the amplitude 
of the leading corrections to scaling with respect to A and D, respectively. 

We also performed simulations of the standard XY model in order to estimate the effect 
of the leading corrections to scaling on the estimates of the critical exponents obtained from 
the FSS analysis. 



D. /?c and the critical value of phenomenological couplings 

In a first step of the analysis we computed R* and the inverse critical temperature /3c at 
A = 2.1 and D = 1.03 respectively. 

For A = 2.1 and D = 1.03 we simulated sc lattices of linear size L from 4 to 16 and 
L = 18, 20, 22, 24, 26, 28, 32, 36, 40, 48, 56, 64, and 80. For all lattice sizes smaller than 
L = 24 we performed 10^ measurements, except for the dd-XY model at L = 16 where 
approximately 1.5 x 10^ measurements were performed. The statistics for the larger lattices 
is given in Table |l[ A measurement was performed after an update cycle as discussed in 

App. [7n. 



Instead of computing R* and Pc from two lattice sizes as discussed in Sec. |II B 2| , we 
perform a fit with the ansatz 

R* = R{L,(3,), (33) 

where R* and /3c are free parameters. We compute R{L, /3) using its third-order Taylor 
expansion 
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A "DT 17 TT 


The number of measurements/ 1000 for the < 


model at A = 2.1 and for the dd-XY 


model at D = 


1.03 for Hnear lattice sizes L > 24. 




L 




dd-XY 


24 


100,000 


112,740 


26 


28,475 


60,900 


28 


101,235 


80,645 


32 


48,005 


55,560 


36 


18,220 


22,880 


40 


21,120 


27,115 


48 


9,735 


16,005 


56 


6,550 


7,495 


64 


4,725 


5,240 


80 


585 


780 



R* = R{L, + d,{L, (3,) {(3, -l3s) + \ d^iL, - (3,)' + ^ ds{L, f3,){f3, - (34) 

/ 

where f3s is the f3 at which the simulation was performed, and R, di, d2, and d^ are determined 
in the MC simulation. 

First, we perform fits for the two models separately. We obtain consistent results for R* 
for all four choices of phenomenological couplings. In order to obtain more precise results for 
(5c and R*, we perform joint fits of both models. Here, we exploit universality by requiring 
that R* takes the same value in both models. Hence, such fits have three free parameters: 
R* and the two values of (3c- In the following we shall only report the results of such joint 
fits. 

Let us discuss in some detail the results for R = Za/Zp that are summarized in Table 
PI . In each fit, we take all data with Lmin ^ L < Lmax into account. For L^ax = 80, the 
X^/d.o.f. becomes approximately one starting from Lmin = 15. However, we should note 
that a x^/d.o.f. close to one does not imply that the systematic errors due to corrections 
that are not taken into account in the ansatz are negligible. 

Our final result is obtained from the fit with Lmin = 28 and Lmax = 80. The systematic 
error is estimated by comparing this result with that obtained using Lmin = 10 and Lmax = 
28. The systematic error on (3c is estimated by the difference of the results from the two fits 
divided by 2.8^'^ — 1, where 2.8 is the scale factor between the two intervals and yt~y-i ~ 2.3. 
Estimating the systematic error by comparison with the interval Lmin = 14 and Lmax = 40 
leads to a similar result. We obtain (3c = 0.5091507(6) [7] for the 0^ model at A = 2.1 
and Pc = 0.5627975(7) [7] for the dd-XY model at D = 1.03. In parentheses we give the 
statistical error and in the brackets the systematic one. Our final result for the critical ratio 
of partition functions is (Za/Zp)* = 0.3202(1) [5]. Here the systematic error is computed by 
dividing the difference of the results of the two fits by 2.8°'^ — 1. 

We repeat this analysis for .^2nd, f^4 and Uq. The final results are summarized in Table 

m 

Next we compute Pc at additional values of A and D. For this purpose we simulated 
lattices of size L = 96 and compute Pc using Eq. (pOf). We use only R = Za/Zp with the 



13 



TABLE III. Joint fits of the Za/Zp data of the (/)^ model and the dd-XY model with the ansatz 
(p^. All lattice sizes Lmin ^ L < Lmax are used in the fit. In column four we give the results of 
the fits for (3c of the cp'^ model at A = 2.1 and in column five the results for (3c of the dd-XY model 
at D = 1.03. Finally, in column six we give the results for the fixed-point value (Za/Zp)*. The 
final results and an estimate of the systematic errors are given in the text. 



r 

-^min 


T 

-'^max 


X /d.o.f. 


(3c, A = 2.1 




Pc, = 1-03 


[Za/Zp) 


11 


80 


3.25 


0.50915354(33) 




0.56280014(35) 


0.319794(25) 


13 


80 


2.48 


0.50915287(35) 




0.56279938(38) 




15 


80 


1.06 


0.50915192(38) 




0.56279834(41) 


0.320019(35) 


20 


80 


0.91 


0.50915142(46) 




0.56279784(49) 


0.320093(52) 


24 


80 


0.89 


0.50915109(53) 




0.56279740(56) 


0.320162(72) 


28 


80 


0.73 


0.50915074(63) 




0.56279747(66) 


0.320195(102) 


32 


80 


0.85 


0.50915065(75) 




0.56279746(82) 


0.320208(149) 


10 


28 


3.47 


0.50915783(53) 




0.56280463(59) 


0.319524(32) 


14 


40 


1.78 


0.50915337(48) 




0.56279981(53) 


0.319877(39) 


TABLE IV. Summary of the final results for Pc and R* . In column one 


the choice of the 


phenomenological coupling R is listed. In columns two and three we report our estimates of (3c for 


the (j)'^ model at A = 


: 2.1 and the dd-XY model at L> = 


: 1.03, 


and finally in column four the result 


for the fixed-point value of the phenomenological coupling. 


Note that the estimates of (3c, based 


on the four different choices of R, are 


consistent within error 


bars. 




R 




/?c, A = 2.1 






= 1.03 


R* 


ZaI Zp 




0.5091507(6) [7] 


0.5627975(7) [7] 


0.3202(1)[5] 


C2nd/ 




0.5091507(7) [3] 


0.5627971(7) [2] 


0.5925(1)[2] 


C4 




0.5091495(9) [10] 


0.5627972(10)[11] 


1.2430(1)[5] 






0.5091498(9) [15] 


0.5627976(10) [15] 


1.7505(3)[25] 


TABLE V. Estimates of f3c from simulations of 96^ 


lattices. (3c is obtained from Eq. (p0|) using 


Za 1 Zp as 


phenomenological coupling. 


In parentheses we give the statistical error and in brackets 


the error 


due to the 


error on [Za/Zp)* 


. "stat" gives (the number of measurements)/1000. 


model 




A ; D 




stat 




Pc 






2.07 




545 




0.5093853(16) [8] 






2.2 




510 




0.5083366(16) [8] 


dd-XY 




0.9 




720 




0.5764582(15) [9] 


dd-XY 




1.02 




1,215 




0.5637972(12) [9] 


dd-XY 




1.2 




665 




0.5470377(17) [9] 



above-reported estimate [Za/Zp)* = 0.3202(6). The results are summarized in Table 0. 
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E. Eliminating leading corrections to scaling 



In this subsection we determine A* and D*. For this purpose, we compute the correction 
amphtude C3 for various choices of Ri and R2 for the 0^ model at A = 2.1 and the dd- 
XY model a.t D = 1.03. In order to convert these results into estimates of A* and D*, 
we determine the derivative of the correction amplitude C3 with respect to A (resp. D) at 
A = 2.1 (resp. D = 1.03). We also simulated the XY model in order to obtain estimates 
of the residual systematic error due to the leading corrections to scaling. Note that, in the 
following, we always use as the value of Ri j in Eq. (0) the estimates of R* given in Table 

m 

1. Derivative of the correction amplitude with respect to X or D 

For this purpose we simulated the dd-XY model at D = 0.9 and D = 1.2 on lattices of 
size L = 5, 6, 7, 8, 9, 10, 12 and 16. The 0^ model was simulated at A = 2.0 and A = 2.2 
on lattices of size L = 3, 4, 5, 6, 7, 8 and 9. In the case of the dd-XY model we performed 
100 X 10^ measurements for each parameter set. In the case of the (p^ model 250 x 10^ 
measurements were performed. 

In the following we discuss only the dd-XY model, since the analysis of the 0"^ data is 
performed analogously. 

In Refs. P^J^ it was observed that subleading corrections to scaling cancel to a large 
extent when one considers the difference of R at close-by values of A. In order to get an idea 
of the size of the corrections, we report in Table ^ 

ARL'-' = {R\d=o.9 - R\d=i.2) L'-' (35) 

for various choices of Ri and i?2- We see that this quantity varies little with L in all cases. In 
the case of -Ri = Za/Zp and R2 = U^, ARL^'^ is already constant within error bars starting 
from L = 5. 

In order to compute C3, see Eq. (^3]), we need ARL'^-^ to be as flat as possible and 
especially large compared to the statistical errors. Looking at Table ^I|, we see that 
the two combinations -Ri = C,2nd/L, R2 = Za/Zp and Ri = f/4, R2 = Uq are unfavorable 
compared with the other four combinations. 

In order to see whether we can predict the exponent 00, we perform a fit with the ansatz 

AR = kL-'^, (36) 

with k and uj as free parameters. From f/4 at Za/Zp = 0.3202 we get u = 0.795(9) with 
X^/d.o.f. = 0.66, using all available data. This value is certainly consistent with field- 
theoretical results. Note however, that we would like to vary the range of the fit in order to 
estimate systematic errors. For this purpose more data at larger values of L are needed. 
In the following we need estimates of 

dD 
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(37) 

D=i.n.s 



TABLE VI. The quantity (-R|d=o.9 - R\d=i.2)L^'^ for the dd-XY model. In the top row we 
give the choice of Ri and R2- For instance, C/4 at {Za/Zp)f means that Ri = Za/Zp and R2 = U^. 

L U4, at {Za/Zp)f Uq at (Zg/Zp)/ Uj at (g2nd/-^)/ Ue at (g2nd/-^)/ Zg/Zp at (6nd/-^)/ ^6 at Ujj 



5 0.0366(2) 


0.1294(7) 


0.0400(2) 


0.1404(7) 


0.0057(2) 


0.0039(1 


6 0.0365(2) 


0.1297(8) 


0.0409(3) 


0.1442(9) 


0.0069(2) 


0.0042(1 


7 0.0368(3) 


0.1312(9) 


0.0415(3) 


0.1469(10) 


0.0073(3) 


0.0046(2 


8 0.0369(3) 


0.1312(10) 


0.0421(3) 


0.1489(11) 


0.0081(3) 


0.0046(2 


9 0.0366(4) 


0.1301(12) 


0.0415(3) 


0.1468(13) 


0.0076(3) 


0.0044(2 


10 0.0368(4) 


0.1311(13) 


0.0419(4) 


0.1483(14) 


0.0078(4) 


0.0045(2 


12 0.0372(4) 


0.1324(15) 


0.0427(5) 


0.1511(17) 


0.0084(4) 


0.0045(3 


16 0.0360(6) 


0.1286(20) 


0.0411(7) 


0.1460(22) 


0.0078(5) 


0.0046(4 



TABLE VII. Estimates for dc3/dD at D = 1.03 (dd-XY) and dcs/dA at A = 2.1 (0^^). In the 
first row we give the combination of Ri and i?2- 



Model 


C/4 at (Za/Zp) f 


Ue at {Za/Zp) f 


Ua at {^2nd/ L)f 


Uq at {^2nd/ L)f 


dd-XY 


-0.122 


-0.435 


-0.140 


-0.495 




-0.0490 


-0.175 


-0.0546 


-0.194 



and of the corresponding quantity for the 0^ model, in order to determine D* and A*. We 
approximated this derivative by a finite difference between D = 0.9 and D = 1.2. The 
coefficient C3 is determined by fixing uj = 0.8. Our final result is the average of the estimates 
for L = 10, 12 and 16 in Table |V1[ In a similar way we proceed in the case of the 0^ model, 
averaging the L = 8,9 results. The results are summarized in Table |VII| . We make no 
attempt to estimate error bars. Sources of error are the finite difference in D, subleading 
corrections, the error on u and the statistical errors. Note however that these errors are 
small enough to be neglected in the following. 



2. Finding R* , A* and D* 

For this purpose we fit our results at D = 1.03 and A = 2.1 with the ansatz 

R = R* + 0^1-'^, (38) 

where we fix = 0.8. We convinced ourselves that setting u = 0.75 or 0.85 changes the 
final results very little compared with statistical errors and errors caused by subleading 
corrections. We perform joint fits, by requiring R* to be equal in both models. 

The results of the fits for four different combinations of Ri and R2 are given in Table 



VIII, where we have already translated the results for C3 into an estimate of A* and D*, by 



usmg 



X* = 2.1--c,[^] (39) 
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TABLE VIII. Results of fits with the ansatz (|3q). The coefficients C3 are converted into D* 
and A* using Eq. (^). All data with Lmin ^ L < L^^,^ are fitted. 







xVd.o.f. 


R* 




A* 


D* 








{Za/Zp)f 


= 0.3202 and R2 


= C/4. 






8 


80 


1.55 




1.24303(2) 




2.077(4) 


1.020(2) 


12 


80 


1.01 




1.24304(4) 




2.071(8) 


1.022(3) 


16 


80 


1.07 




1.24308(6) 




2.057(14) 


1.019(5) 


20 


80 


1.12 




1.24301(8) 




2.073(22) 


1.028(9) 


8 


40 


1.62 




1.24304(3) 




2.077(4) 


1.020(2) 


10 


40 


1.02 




1.24305(3) 




2.070(6) 


1.020(2) 






^1,/ = 


{Za/Zp)f 


= 0.3202 and R2 








8 


80 


2.15 




1.75156(8) 




2.006(4) 


0.990(2) 


12 


80 


1.15 




1.75126(13) 




2.018(7) 


1.000(3) 


16 


80 


1.22 




1.75120(19) 




2.017(13) 


1.003(5) 


20 


80 


1.19 




1.75085(27) 




2.043(21) 


1.015(8) 


8 


40 


2.21 




1.75160(9) 




2.004(4) 


0.989(2) 


10 


40 


1.24 




1.75143(11) 




2.010(6) 


0.994(2) 






^1,/ = 


(?2nd/ -^)/ 


= 0.5925 and R2 








8 


80 


4.01 




1.24352(3) 




1.977(4) 


0.987(2) 


12 


80 


1.19 




1.24322(4) 




2.031(8) 


1.010(3) 


16 


80 


1.29 




1.24314(6) 




2.049(14) 


1.019(5) 


20 


80 


1.13 




1.24299(9) 




2.083(23) 


1.035(9) 


8 


40 


4.19 




1.24355(3) 




1.973(4) 


0.985(2) 


10 


40 


1.49 




1.24335(4) 




2.006(6) 


1.000(2) 






^1,/ = 


iC2nd /L)f 


— n ^091^ anA 

— v.Oij^o anu. 1x2 


— TTr- 






8 


80 


6.62 




1.75323(10) 




1.915(4) 


0.961(2) 


12 


80 


1.55 




1.75189(14) 




1.985(7) 


0.991(3) 


16 


80 


1.51 




1.75142(22) 




2.013(13) 


1.004(5) 


20 


80 


1.24 




1.75078(32) 




2.055(22) 


1.024(8) 


8 


40 


7.02 




1.75336(10) 




1.911(4) 


0.959(2) 


10 


40 


2.21 




1.75248(12) 




1.953(6) 


0.978(2) 



for the 0^ model and the analogous formula for the dd-XY model, and the results of Table 



A x^/d.o.f. close to 1 is reached for L^i^ = 10 and Lmax = 80 in the case of f/4 at 
{Za/Zp)f = 0.3202. This has to be compared with Ljam = H, H, and 14 in the case of Uq 
at {Za/Zp)f = 0.3202, f/4 at (6nd/^)/ = 0.5925 and Uq at (6nd/^)/ = 0.5925. 

This indicates that U4 at {Za/Zp)f = 0.3202 has the least bias due to subleading correc- 
tions to scaling. Therefore we take as our final result A* = 2.07 and D* = 1.02 which is the 
result of Lmin = 12 and Lmax = 80 in Table [VIII| . Starting from Lmin = 20 all results for A* 
and D* are within 2a of our final result quoted above. 

Our final results are A* = 2.07(5) and D* = 1.02(3). The error bars are such to include 
all results in Table |Vlli| with Lmin = 20 and L^a.x = 80, including the statistical error, and 
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TABLE IX. Corrections to scaling in the standard XY model for R with R2 = f/4 and 
R^ j = {Za/Zp)f = 0.3202. We use the ansatz (||) with uj = 0.8 fixed. 



X /d.O.l R* C3 

12 64 1.78 1.2432(1) -0.1120(7) 

16 64 0.73 1.2430(1) -0.1087(13) 

20 64 0.38 1.2427(2) -0.1048(22) 

24 64 0.24 1.2429(2) -0.1083(34) 

12 32 2.32 1.2433(1) -0.1124(8) 



therefore should take into proper account all systematic errors. 

From these results, it is also possible to obtain a conservative upper bound on the coeffi- 
cient C3 for A = 2.1 and D = 1.03. Indeed, using the estimates of A* and D* and their errors, 
we can obtain the upper bounds |2.1 - A*| < AA = 0.08 and |1.03 - D*\ < AD = 0.04. 
Then, we can estimate |c3(A = 2.1)| < AX^dcs/dX), and analogously Ic^^D = 1.03)| < 
AD{dc3/dD). For at {Za/Zp)f = 0.3202, using the results of Table we have 

|C3(A = 2.1)1 < 0.004, |c3(L> = 1.03)1 < 0.005. (40) 



3. Corrections to scaling in the standard XY model 

We simulated the standard XY model on lattices with linear sizes L = 6, 8, 10, 12, 16, 
18, 20, 22, 24, 28, 32, 48, 56, and 64 at f3, = 0.454165, which is the estimate of of Ref. H]. 



Here, we used only the wall-cluster algorithm for the update. In one cycle we performed 12 
wall-cluster updates. For L < 16 we performed 10® cycles. For lattice sizes 16 < L < 64, we 
spent roughly the same amount of CPU time for each lattice size. For L = 64 the statistics 
is 2.73 X 10^ measurements. 

We determine the amplitude of the corrections to scaling for R with Ri j = {Za/Zp)f = 
0.3202 and R2 = U4. Other choices lead to similar results. We fit our numerical results with 



the ansatz (p8D, where we fix = 0.8. The results are given in Table IX. Note that the 



results for R* are consistent with the result obtained from the joint fit of the two improved 



models. In Table ^TTl| we obtained, e.g., R* = 1.24301(8) with L^in = 20 and L^ax = 80. 

Corrections to scaling are clearly visible, see Fig. |l|. From the fit with Lmin = 20 and 
-^max = 64 we obtain C3 = —0.1048(22). For the following discussion no estimate of the 
possible systematic errors of C3 is needed. Comparing with Eq. (^Of ), we see that in the 
(approximately) improved models the amplitude of the leading correction to scaling is at 
least reduced by a factor of 20. Note, that even if this result was obtained by considering 
a specific observable, [4 at fixed Za/Zp, the universality of the ratios of the subleading 
corrections implies the same reduction for any quantity. In the following section we will use 
this result to estimate the systematic error on our results for the critical exponents. 
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FIG. 1. Corrections to scaling for the dd-XY model at D = 0.9, 1.03, and 1.2, and for the 
standard XY model. We plot R with Ri j = {Za/Zp)f = 0.3202 and R2 = {74 as a function of the 
lattice size. The dotted lines should only guide the eye. 

F. Critical exponents from finite-size scaling 



As discussed in Sec. |IIB| , we may use the derivative of phenomenological couplings taken 
at f3f in order to determine yt- Given the four phenomenological couplings that we have 
implemented, this amounts to 16 possible combinations. In the following we will restrict the 
discussion to two choices: in both cases we fix Pf by {Za/Zp)f = 0.3202. At /?/ we consider 
the derivative of the Binder cumulant and the derivative of Za/Zp. In Table ^ we summarize 
the results of the fits with the ansatz 



dR 
'dp 



(41) 



for the (f)^ model at A = 2.1, the dd-XY model sX D = 1.03, and the standard XY model. 

We see that for the same Lmin and Lmax the statistical error on the estimate of v obtained 
from the derivative of Za/Zp is smaller than that obtained from the derivative of f/4. On the 
other hand, for the two improved models, scaling corrections seem to be larger for Za/Zp 
than for f/4. 

In the case of Za/Zp, for both improved models, the result of the fit for v is increasing 
with increasing Lmin. In the case of the Binder cumulant, it is increasing with Lmin for the 0^ 
model and decreasing for the dd-XY model. The fact that scaling corrections affect the two 
quantities and the two improved models in a quite different way suggests that systematic 
errors in the estimate of v can be estimated from the variation of the fits presented in Table 
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TABLE X. Fit results for the critical exponent v obtained from the ansatz (^). In all cases /9/ 
is fixed by Za/Zpj = 0.3202. We analyze the (j)^ model at A = 2.1, the dd-XY model at D = 1.03, 



and the standard XY model. We consider the slope of the Binder cumulant C/4 and of the ratio of 
partition functions Za/Zp. We included all data with L^\a < L < Lmax into the fit. 







xVd.o.f. 


V 


(j)^ model: derivative of U4 


7 


80 


1.17 


0.67168(12) 


9 


80 


0.79 


0.67188(15) 


11 


80 


0.85 


0.67181(19) 


16 


80 


0.98 


0.67192(34) 


(f)'^ model: derivative of Za/Zp 


12 


80 


3.01 


0.67042(9) 


16 


80 


1.61 


0.67104(15) 


20 


80 


1.04 


0.67139(22) 


24 


80 


0.54 


0.67194(32) 


dd-XY model: derivative of C/4 


7 


80 


2.06 


0.67258(12) 


9 


80 


1.13 


0.67216(15) 


11 


80 


1.19 


0.67209(19) 


16 


80 


0.97 


0.67154(31) 


dd-XY model: derivative of Za/Zp 


12 


80 


1.89 


0.67017(9) 


16 


80 


1.60 


0.67046(14) 


20 


80 


0.79 


0.67099(21) 


24 


80 


0.80 


0.67113(30) 


XY model: derivative of C/4 


12 


64 


4.48 


0.66450(28) 


16 


64 


1.30 


0.66618(42) 


20 


64 


0.54 


0.66740(63) 


XY model: derivative of Za/Zp 


12 


64 


1.33 


0.67263(13) 


16 


64 


0.69 


0.67300(19) 


20 


64 


0.30 


0.67325(30) 


24 


64 


0.25 


0.67327(41) 



As our final result we quote v = 0.6716(5) which is consistent with the two results from 
Za/Zp at Ljnin = 24 and with the results from U4 at L^in = 16. 

In the case of the standard XY model, the derivative of f/4 requires a much larger Lmin 
to reach a small x^/d.o.f. than for the improved models. For the derivative of Za/Zp instead 
a x^/d.o.f. fti 1 is obtained for an Lmin similar to that of the improved models. Note that 
the result for u from the derivative of U4 for Lmin = 16 is by several standard deviations 
smaller than our final result from the improved models, while the result from the derivative 
of Za/Zp is by several standard deviations larger! Again we have a nice example that a 
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TABLE XI. Results for the critical exponent rj from the FSS of the maj 


;netic susceptibility. 


Fits with ansatz ( 


42). All data with Lmin 


< L < L^ioi are taken into account. 




-^min 




XV d.o.f. 


V 


(/)^ model 


20 


80 


2.44 


0.0371(1) 


24 


80 


0.73 


0.0375(1) 


28 


80 


0.94 


0.0375(2) 


32 


80 


0.41 


0.0378(3) 


dd-XY model 


20 


80 


1.88 


0.0371(1) 


24 


80 


1.19 


0.0373(1) 


28 


80 


1.52 


0.0374(2) 


32 


80 


1.24 


0.0376(2) 


XY model 


20 


64 


7.92 


0.0325(2) 


24 


64 


1.81 


0.0344(2) 


28 


64 


0.27 


0.0340(3) 


32 


64 


0.06 


0.0342(4) 



X^/d.o.f. ~ 1 does not imply that systematic errors due to corrections that have not been 
taken into account in the fit are small. 

Remember that in improved models the leading corrections to scaling are suppressed at 
least by a factor of 20 with respect to the standard XY model. Since the range of lattice sizes 
is roughly the same for the XY model and for the improved models, we can just divide the 
deviation of the XY results from v = 0.6716(5) by 20 to obtain an estimate of the possible 
systematic error due to the residual leading corrections to scaling. For the derivative of 
Za/Zp we end up with 0.0001 and for the derivative of U4 with 0.0003. 

We think that these errors are already taken into account by the spread of the results for 
u from the derivatives of 1/4^ and Za/Zp and the two improved models. Therefore, we keep 
our estimate u = 0.6716(5) with its previous error bar. 

Next we compute the exponent rj. For this purpose we study the finite-size behavior of 
the magnetic susceptibility at Pj. In the following we fix Pf by Rij = (Za/Zp) f = 0.3202. 
Other choices for Ri j give similar results. 

In a first attempt we fit the data of the two improved models and the standard XY model 
to the simple ansatz 

X\,, = CL'-^. (42) 



The results are summarized in Table XI 



For all three models rather large values of Lmin are needed in order to reach a x /d.o.f. 
close to one. In all cases the estimate of 77 is increasing with increasing Lmin- For Lmin = 24 
the result for rj from the standard XY model is lower than that of the improved models by an 
amount of approximately 0.0030 . Therefore, the systematic error due to leading corrections 
on the results obtained in the improved models should be smaller than 0.0030/20 = 0.00015. 
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TABLE XII. Results for the critical exponent rj from the FSS of the magnetic susceptibility. 
Fits with ansatz (|4^). All data with Lmin ^ L < Lmax are taken into account. 







fit with K = 0.0 


fit with K 


= 0.2 






XV d.o.f. 


r] 


XV d.o.f. 




(j)'^ model 


8 


80 


0.72 


0.0386(1) 


1.16 


0.0391(1) 


10 


80 


0.68 


0.0385(1) 


1.27 


0.0388(1) 


12 


80 


0.75 


0.0385(1) 


0.81 


0.0388(1) 


14 


80 


0.84 


0.0386(2) 


0.92 


0.0388(2) 


16 


80 


0.72 


0.0384(2) 


0.73 


0.0386(2) 


20 


80 


0.88 


0.0384(3) 


0.88 


0.0385(4) 


dd-XY model 


8 


80 


1.85 


0.0387(1) 


3.06 


0.0391(1) 


10 


80 


0.95 


0.0384(1) 


1.15 


0.0388(1) 


12 


80 


0.99 


0.0384(1) 


1.04 


0.0386(1) 


14 


80 


0.94 


0.0384(2) 


1.03 


0.0386(2) 


16 


80 


0.85 


0.0383(2) 


1.14 


0.0384(2) 


20 


80 


0.90 


0.0381(4) 


0.90 


0.0382(4) 


XY model 


12 


64 


0.64 


0.0350(2) 






16 


64 


0.48 


0.0353(3) 






20 


64 


0.37 


0.0358(5) 







Given this tiny effect, it seems plausible that, for the improved models, the increase of the 
estimate of r/ with increasing Lmin is caused by subleading corrections. Therefore, we consider 
0.0375, which is the result of the fit with Lmin = 24 in the 0^ model, as a lower bound of rj. 

Finally, we perform a fit which takes into account the analytic background of the magnetic 
susceptibility. In Ref. P|, it was shown that the addition of a constant term to Eq. (^21) 
leads to a small x^/d.o.f. already for small Lmin < 10. Similar results have been found for 
the Ising universality class. This ansatz is not completely correct, since it does not take into 
account corrections proportional to with y —1.8, which formally are more important 
than the analytic background. However, the difference between these exponents is small, 
and a four-parameter fit is problematic Therefore, we decided to fit our data with the ansatz 

x\p^=cL^~^ + hL\ (43) 

with K fixed to 0.0 and to 0.2. The difference between the results of the fits with the two 
values of k will give an estimate of the systematic error of the procedure. Results for all 
three models are summarized in Table PCII] . 

The value x^/d-o.f. is close to one for Lmin = 8 for the (p'^ model and Lmin = 10 for 
the dd-XY model, and it does not allow to discriminate between the two choices of k. The 
values of r] are rather stable as Lmin is varied, although there is a slight trend towards smaller 
results as Lmin increases; the trend seems to be stronger for /t = 0.2. Moreover, the results 
from the two models are in good agreement. 
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The fits for the XY model also give a good x^/d-o.f. for Lmin > 12; the value of 77 is 
however much too small, and shows an increasing trend. We can estimate from the difference 
between the XY model and the improved models at Lmin = 16 that the error on the value of 
rj obtained from improved models, induced by residual leading scaling corrections, is smaller 
than 0.003/20 = 0.00015. 



From the results for the improved models reported in Table PCII] , one would be tempted 
to take 1] = 0.0384 as the final result. However, as we can see from the results for the XY 
model, we should not trust blindly the good x^/d.o.f. of these fits. Taking into account 
the decreasing trend of the values of 1] for the improved models, we assign the conservative 
upper bound r] < 0.0385. By combining it with the lower bound obtained from ansatz (^2[), 
we obtain our final result 

0.0375 <ri< 0.0385, i.e., r] = 0.0380(5). (44) 



III. HIGH-TEMPERATURE DETERMINATION OF CRITICAL EXPONENTS 

In this Section we report the results of our analyses of the HT series. The details are 
reported in App. 0. 

We compute 7 and u from the analysis of the HT expansion to Olf]"^^) of the magnetic 
susceptibility and of the second-moment correlation length. In App. |B 2| we report the details 
and many intermediate results so that the reader can judge the quality of our results without 
the need of performing his own analysis. This should give an idea of the reliability of our 
estimates and of the meaning of the errors we quote, which depend on many somewhat 
arbitrary choices and are therefore partially subjective. 

We analyze the HT series by means of integral approximants (lA's) of first, second, 
and third order. The most precise results are obtained biasing the value of (3c, using its 
MC estimate. We consider several sets of biased lA's, and for each of them we obtained 
estimates of the critical exponents. These results are reported in App. [B 2| . All sets of lA's 
give substantially consistent results. Moreover, the results are also stable with respect to the 
number of terms of the series, so that there is no need to perform problematic extrapolations 
in the number of terms in order to obtain the final estimates. The error due to the uncertainty 
on A* and D* is estimated by considering the variation of the results when changing the 
values of A and D. 

Using the intermediate results reported in App. P 2| , we obtain the estimates of 7 and u 
shown in Table pClll[ We report on 7 and u three contributions to the error. The number 
within parentheses is basically the spread of the approximants at the central estimate of A* 
{D*) using the central value of (3c- The number within brackets is related to the uncertainty 
on the value of (3c and is estimated by varying (3c within one error bar at A = A* or D = D* 
fixed. The number within braces is related to the uncertainty on A* or D*, and is obtained by 
computing the variation of the estimates when A* or D* vary within one error bar, changing 
correspondingly the values of (3c- The sum of these three numbers should be a conservative 
estimate of the total error. 
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TABLE XIII. Estimates of the critical exponents obtained from the analysis of the HT expan- 
sion of the improved 0^ lattice Hamiltonian and dd-XY model. 



7 


V 




a 


(j)'^ Hamiltonian 1.31780(10)[27]{15} 
dd-XY model 1.31748(20)[22]{18} 


0.67161(5)[12]{10} 
0.67145(10)[10]{15} 


0.0380(3){1} 
0.0380(6){2} 


-0.0148(8) 
-0.0144(10) 



We determine our final estimates by combining the results for the two improved Hamil- 
tonians: we take the weighted average of the two results, with an uncertainty given by the 
smallest of the two errors. We obtain for 7 and v 

7 = 1.3177(5), (45) 
V = 0.67155(27), (46) 

and by the hyperscaling relation a = 2 — Su 

a = -0.0146(8). (47) 

Consistent results, although significantly less precise (approximately by a factor of two), are 
obtained from the IHT analysis without biasing jSc (see App. p 2| ). 

From the results for 7 and z/, we can obtain rj by the scaling relation 7 = {2 — rfjv. This 
gives rj = 0.0379(10), where the error is estimated by considering the errors on 7 and u as 
independent, which is of course not true. We can obtain an estimate of 77 with a smaller, 
yet reliable, error by applying the so-called critical-point renormalization method (CPRM) 
(see, e.g., Refs. and references therein) to the series of x and The results are reported 
in Table pClll|. We report two contributions to the error on 77, as discussed for 7 and z/; the 



uncertainty on Pc does not contribute in this case. Our final estimate is 

r] = 0.0380(4). (48) 
Moreover, using the scaling relations we obtain 

5=^= 4.780(2), (49) 
/3 = ^(l + r/) = 0.3485(2), (50) 
where the error on (3 has been estimated by considering the errors of u and rj as independent. 

IV. THE CRITICAL EQUATION OF STATE 
A. General properties of the critical equation of state of XY models 

We begin by introducing the Gibbs free-energy density 

G{H) = ^\ogZ{H), (51) 
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and the related Helmholtz free-energy density 



J^{M) = M ■ H -G{H), (52) 

— * — * 

where V is the volume, M the magnetization density, H the magnetic field, and the depen- 
dence on the temperature is understood in the notation. A strictly related quantity is the 
equation of state which relates the magnetization to the external field and the temperature: 

(53) 

dM 

In the critical limit, the Helmholtz free energy obeys a general scaling law. Indeed, for t ^ 0, 
\M\ 0, and t|M|-V/3 fixed, it can be written as 

AJ^ = J^(M) - J^regiM) ~ t2-°^(|M|t-^), (54) 

where Tj-egiM) is a regular background contribution. The function JF is universal apart from 
trivial rescalings. 

The Helmholtz free energy is analytic outside the critical point and the coexistence curve 
(Griffiths' analyticity). Therefore, it has a regular expansion in powers of \M\ for t > fixed, 
which we write in the form 

AT = \m\' + J2 ^"''j^92,^"'^ (55) 

where m = 1/^, ^ is the second- moment correlation length, is a renormalized magnetiza- 
tion, and g2j are the zero-momentum 2j-point couplings. By performing a further rescaling 
if = Wb^l'^zj ^fgl^ the free energy can be written as P7[ 

Tfi 

AT=—A{z), (56) 

where 

^W-i^^ + ^^' + E^^r,,.^'. (57) 

Note that z oc \M\t~^ for t — 0, so that Eq. ( pTj ) is nothing but the expansion of T{\M\t~^^ 
for |M|t~^ 0. Correspondingly, by using the scaling relation /3(1 + 5) = 2 — a, we obtain 
for the equation of state 

E oc t^^F{z), (58) 

with 
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Because of Griffiths' analyticity, J^{M) has also a regular expansion in powers of t for \M\ 
fixed. Therefore, 

oo oo 
fc=0 fc=0 

where, because of Eq. (H), the coefficients Tk{M) scale as |M|-('=+°-2)//3 for |M| ^ 0. From 
this expression we immediately obtain the large- 2; expansion of -^(-2), 

F{z) = z^^F^z-^'^, (61) 

fc=0 

where we have used again (3{l + 5) = 2 — a. The function F{z) is defined only for t > 0, and 
thus, in order to describe the low-temperature region t < 0, one should perform an analytic 
continuation in the complex t plane [^,0]. The coexistence curve corresponds to a complex 
Zq = \zQ\e~'^'^^ such that -^(2:0) = 0. Therefore, the behavior near the coexistence curve is 
related to the behavior of F{z) in the neighborhood of zq. The constants and |zo| can 
be expressed in terms of universal amplitude ratios, by using the asymptotic behavior of the 
magnetization along the critical isotherm and at the coexistence curve: 

(r<+\{35-l)/2 

rpoo _ \^ > /gn\ 

^0 (5C"=)5(-c+)('5-i)/2' y^^' 

\z,\' = Rt ^ -CtB'/{C^f, (63) 

where the critical amplitudes are defined in App. |C|. 

The function F{z) provides in principle the full equation of state. However, it has the 
shortcoming that temperatures t < correspond to imaginary values of the argument. It is 
thus more convenient to use a variable proportional to t\M\~^^^, which is real for all values 
of t. Therefore, it is convenient to rewrite the equation of state in a different form, 

H = M\M\^-^f{x), X (xt\M\-'^^^, (64) 

where f{x) is a universal scaling function normalized in such a way that /(—I) = and 
/(O) = 1. The two functions f{x) and F{z) are clearly related: 



z 



-'F{z)=F^f{x), z=\zo\x-^. (65) 



It is easy to reexpress the results we have obtained for F{z) in terms of x. The regularity 
of F[z) for 2; — > implies a large-x expansion of the form 



00 



/(a;)=a;^5^/ra;-^"^ (66) 



n=0 



The coefficients can be expressed in terms of r2„ using Eq. (E^. We have 



Fo-(2n + l)!' ^^^^ 
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where we set r2 = r4 = 1. In particular, using Eqs. (|62D and (|6^), 

= R-^\ (68) 

where is defined in App. |C[ Finally, we notice that Griffiths' analyticity implies that 
f{x) is regular for x > —1. In particular, it has a regular expansion in powers of x. The 
corresponding coefficients are easily related to those appearing in Eq. (|6TD. 

We want now to discuss the behavior of f{x) for x —1, i.e., at the coexistence curve. 
General arguments predict that at the coexistence curve the transverse and longitudinal 
magnetic susceptibilities behave respectively as 

M dM 
XT = ^^=dH^ ^ ■ ^^^^ 

In particular the singularity of xl for t < and H ^ is governed by the zero-temperature 
infrared-stable Gaussian fixed point |^9|-^, leading to the prediction 

f{x) ~ c/ (1 + xf for X -1. (70) 



The nature of the corrections to the behavior ([7DD is less clear. It has been conjectured 



5^ , |53|j51| , using essentially e-expansion arguments, that, for y ^ 0, i.e., near the coexistence 
curve, V = 1 + X has a double expansion in powers oi y = HM~^ and y^'^^'^y^. This implies 
that in three dimensions f{x) can be expanded in powers of v at the coexistence curve. On 
the other hand, an explicit calculation to next-to-leading order in the expansion 
shows the presence of logarithms in the asymptotic expansion of f{x) for x — >■ — 1. However, 
they are suppressed by an additional factor of v^logv compared to the leading behavior 

It should be noted that for the A transition in ^He the order parameter is related to 
the complex quantum amplitude of helium atoms. Therefore, the "magnetic" field is not 
experimentally accessible, and the function appearing in Eq. (|6^) cannot be measured di- 
rectly in experiments. The physically interesting quantities are universal amplitude ratios of 
quantities formally defined at zero external field, such as Uq = j A~ ^ for which accurate 
experimental estimates have been obtained. On the other hand, the scaling equation of state 
is physically relevant for planar ferromagnetic systems. 



B. Small- A/ expansion of the equation of state in the high-temperature phase 



Using HT methods, it is possible to compute the first coefficients and appearing 
in the expansion of the Helmholtz free energy and of the equation of state, see Eqs. and 
(^). Indeed, these quantities can be expressed in terms of zero-momentum 2j-correlation 
functions and of the correlation length. 

Details of the analysis of the HT series of (74, rg, rg, and rio are reported in App. [B 3|. 



We obtained the results shown in Table XIV. In Table XV we report our final estimates 



(denoted by IHT), obtained by combining the results of the two models; we also compare 
them with the estimates obtained using other approaches. Note that our final estimate of 
(?4 is slightly larger than the result reported in Ref. ||2^ (see Table pCV| ). The difference is 
essentially due to the different analysis employed here, which should be more reliable. This 
point is further discussed in App. 
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TABLE XIV. Estimates of (74, rg, rg, and rio obtained from the analysis of the HT series for 



the two improved Hamiltonians. Final results will be reported in Table XV 



94 re rs rio_ 

(p^ Hamiltonian 21.15(6) 1.955(20) 1.37(15) -13(7) 

dd-XY model 21.13(7) 1.948(15) 1.47(10) -11(14) 



TABLE XV. Estimates of (74, rg, rg, and rio obtained using the following methods: anal- 
yses of improved HT expansions (IHT), of HT expansions for the standard XY model (HT), of 
fixed-dimension perturbative expansions (d = 3 g-exp.), and of e expansions (e-exp.). A more 
precise determination of rio will be reported in Table XVI| . 





IHT 


HT 


d = 3 g-exp. 


e-exp. 


94 


21.14(6) 
21.05(6) m 


21.28(9) m 
21.34(17) 


21.16(5) i 

21.11 g 


21.5(4) (HJll 


re 


1.950(15) 
1.951(14) ^ 


2.2(6) 11 


1.967 13 


1.969(12) l55,58| 


rs 


1.44(10) 

1.36(9) mi 




1.641 [57| 


2.1(9) |5,58| 


no 


-13(7) 









C. Parametric representations of the equation of state 

In order to obtain a representation of the equation of state that is vaUd in the whole 
critical region, we need to extend analytically the expansion (|59D to the low-temperature 
region t < 0. For this purpose, we use parametric representations that implement the 
expected scaling and analytic properties. They can be obtained by writing [p^-|6T|| 

M = moR'^m{9), 

t = R{i-e^), 

H = hoR'^^hie), (71) 

where ho and mo are normalization constants. The variable R is nonnegative and measures 
the distance from the critical point in the {t, H) plane, while the variable 9 parametrizes the 
displacement along the lines of constant R. The functions m{9) and h{9) are odd and regular 
at ^ = and at ^ = 1. The constants mo and can be chosen so that m{9) = 9 + 0{9^) 
and h{9) = 9 + 0{9^). The smallest positive zero of h{9), which should satisfy ^0 > 1, 
corresponds to the coexistence curve, i.e., to T < Tc and H ^ 0. The singular part of the 
free energy is then given by 

AJ^ = homoR^-''g{9), (72) 
where g{9) is the solution of the first-order differential equation 

(1 - 9^)g'{9) + 2(2 - a)9g{9) = [{1 - 9^)m'{9) + 2f39m{9)] h{9) (73) 
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that is regular at 6 = 1. In particular, the ratio A'^ /A of the specific-heat amplitudes in 
the two phases can be derived by using the relation 



AVA- = (6^ - 1)^-"^. (74) 

The parametric representation satisfies the requirements of regularity of the equation of 
state. Singularities can appear only at the coexistence curve (due, for example, to the 



logarithms discussed in Ref. |54]), i.e., for 6 = 6q. Notice that the mapping (71) is not 



invertible when its Jacobian vanishes, which occurs when 

Y{0) = {1- e^)m\e) + 2f5em{e) = o. (75) 

Thus, parametric representations based on the mapping (^) are acceptable only if Oq < 9i 

where di is the smallest positive zero of the function Y{6). One may easily verify that the 
asymptotic behavior (^) is reproduced simply by requiring that 

h{e) ~ (^0 - for 9 60. (76) 

The functions m{6) and h{6) are related to F{z) by 

z = pm{e){i-ey'' , (77) 

Fizie)) = p{i-e')-''hie), (78) 



where p is a free parameter |^,|16|. In the exact parametric equation the value of p may 
be chosen arbitrarily but, as we shall see, when adopting an approximation procedure the 
dependence on p is not eliminated. In our approximation scheme we will fix p to ensure 
the presence of the Goldstone singularities at the coexistence curve, i.e., the asymptotic 
behavior ([76|). Since z = p6 -\- 0{6^), expanding m{6) and h{6) in (odd) powers of 6, 



n(^) = + ^m2n+i^'"+\ 

n=l 

hi0) = 9 + J2h2n+i0'''^\ (79) 



n=l 



and using Eqs. (|77D and ([78D, one can find the relations among p, m2n+i, ^2n+i and the 
coefficients r2„ of the expansion of F{z). 

Following Ref. p6[, we construct approximate polynomial parametric representations 



that have the expected singular behavior at the coexistence curve pD|- pT] , |5i| (Goldstone 



singularity) and match the known small-z expansion (^). We will not repeat here in full 



the discussion of Ref. ||26|, which should be consulted for more details. We consider two 
distinct schemes of approximation. In the first one, which we denote by (A), h{6) is a 
polynomial of fifth order with a double zero at ^o, and m{6) a polynomial of order (1 + 2n): 

(n ^ 
1 + CiO^' 
i=l J 

h{e) = d{i-e^/elf . (so) 



29 



In the second scheme, denoted by (B), we set 



scheme (B) : m{6) = 6, 

h{9) =9{1- 9'/ 9',)' + I • (81) 



i=l 



Here h{9) is a polynomial of order 5 + 2n with a double zero at ^o- Note that for scheme 
(B) 

Y{9) = 1-9^ + 2(39^, (82) 

independently of n, so that = (1 — 2f3)~^. Concerning scheme (A), we note that the 
analyticity of the thermodynamic quantities for \9\ < 9q requires the polynomial function 
Y{9) not to have complex zeroes closer to the origin than 9o- 

In both schemes the parameter p is fixed by the requirement (^), while 9o and the n 
coefficients Cj are determined by matching the small- 2; expansion of F{z). This means that, 
for both schemes, in order to fix the n coefficients Cj we need to know n + 1 values of r2j, 
i.e., rg, ...r6+2n- As input parameters for our analysis we consider the estimates obtained in 
this paper, i.e., a = -0.0146(8), r] = 0.0380(4), rg = 1.950(15), rg = 1.44(10), = -13(7). 

Before presenting our results, we mention that the equation of state has been recently 
studied by MC simulations of the standard XY model, obtaining a fairly accurate determi- 
nation of the scaling function f{x) [Q. In particular we mention the precise result obtained 
for the universal amplitude ratio (see App. |C|for its definition), i.e., = 1.356(4), and 
for the constant Cf, i.e., Cf = 2.85(7), where Cf is defined in Eq. ([7D|). In the following we 
will take into account these results to find the best parametrization within our schemes (A) 
and (B). 

By using the few known coefficients — essentially rg and rg because the estimate of rio 
is not very precise — one obtains reasonably precise approximations of the scaling function 
F{z) for all positive values of z, i.e., for the whole HT phase up to t = 0. In Fig. ^ we 
show the curves obtained in schemes (A) and (B) with n = 1 that use the coefficients rg 
and rg. The two approximations of F{z) are practically indistinguishable. This fact is not 
trivial since the small- 2; expansion has a finite convergence radius which is expected to be 
\zq\ = (i?^)^/^ ^ 2.7. Therefore, the determination of F{z) on the whole positive real axis 
from its small- 2; expansion requires an analytic continuation, which turns out to be effectively 
performed by the approximate parametric representations we have considered. 

Larger differences between the approximations given by schemes (A) and (B) for n = 1 
appear in the scaling function f{x), which is shown in Fig. ^, especially in the region a; < 0, 
which corresponds to t < and z imaginary. Note that the sizeable differences for x > are 
essentially caused by the normalization of f{x), which is performed at the coexistence curve 
X = —1 and at the critical point x = 0, by requiring /(—I) = and /(O) = 1. Although the 
large-x region corresponds to small values of z, the difference between the two approximate 
schemes does not decrease in the large-x limit due to their slightly different estimates of Ry_ 
(see Table PCV1| ). Indeed, for large values of x 



fix) ~ R-'x\ (83) 
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FIG. 2. The scaling function F{z) 
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FIG. 3. The scaHng function f{x). The MC curve is taken from Ref. p 



In Fig. y we also plot the curve obtained in Ref. |^ by fitting the MC data. 

In Table [XVI| we report the results for some universal ratios of amplitudes. The no- 
tations are explained in App. |C[ The reported errors refer only to the uncertainty of the 
input parameters and do not include the systematic error of the procedure, which may be 
determined by comparing the results of the various approximations. Comparing the results 
for and c/ with the MC estimates of Ref. |6^, we observe that the parametrization (A) 
is in better agreement with the numerical data. This is also evident from Fig. ^. 
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TABLE XVI. Universal ratios of amplitudes computed using a = -0.0146(8), r/ = 0.0380(4), 
rg = 1.950(15) and rg = 1.44(10), rio = —13(7). The numbers of the first four lines correspond to 
central values of the input parameters. The errors reported are only related to the uncertainty on 
the input parameters. Numbers marked with an asterisk are inputs, not predictions. 

[(A) n = 1; r6,r8] [(B) n = 1; re, rg] [(A) n = 2; re, rg, fi^] [(B) n = 2; re, rg, rip] 



p 

f 


2.22974 


2.06825 


2.23092 


2.04 


U 


3.88383 


2.93941 


3.88686 


2.70 


Cl 


-0.0260296 


0.0758028 


-0.0265055 


0.11 


C2 








0.0002163 


0.01 


A+/A- 


1.062(4) 


1.064(4) 


1.062(3) 


1.062(5) 


Rt 


0.355(3) 


0.350(1) 


0.355(2) 


0.354(5) 


Rc 


0.127(6) 


0.115(2) 


0.126(2) 


0.119(8) 


R^ 


1.35(7) 


1.50(2) 


*1.356 


1.45(8) 


R^ 


7.5(2) 


7.92(8) 


7.49(6) 


7.8(3) 


7700 


0.0302(3) 


0.0300(2) 


0.0302(2) 


0.0302(4) 


no 


-10(1) 


-11.9(1.4) 


-10(1) 


*-13(7) 


c/ 


4(2) 


52(20) 


4(2) 





We also consider both schemes with n = 2. If we use rio to determine the next coefficient 
C2, scheme (A) is not particularly useful because it is very sensitive to rio, whose estimate 
has a relatively large error [|26|. This fact was already observed in Ref. |2^, and explained 



by considerations on the more complicated analytic structure. One may instead determine 
C2 by using the MC result = 1.356(4). The estimates of the universal amplitude ratios 
obtained in this way are presented in Table PCVI| . They are very close to the n = 1 case. 



providing additional support to our estimates and error bars. Scheme (B) is less sensitive to 
rio and provides reasonable results if we use rio to fix the coefficient C2 in h{6) and impose 
the consistency condition 9q < 9i. The results are shown in Table pCVl| , where one observes 
that they get closer to the estimates obtained by using scheme (A). 

As already mentioned, the most interesting quantity is the specific-heat amplitude ratio 
/ A~ , because its estimate can be compared with experimental result. Our results for 
A'^ / A^ are quite stable and insensitive to the different approximations of the equation of 
state we have considered, essentially because they are obtained from the function g{9), which 
is not very sensitive to the local behavior of the equation of state, see Eq. (^3]). From Table 
XVll we obtain the estimate 



A-^/A- = 1.062(4). (84) 
In Table pCVll| we compare our result (denoted by IHT-PR) with other available es- 



timates. Note that there is a marginal disagreement with the result of Ref. [^, i.e 



y4+/A^ = 1.055(3), which was obtained using the same method but with different input 
parameters: a = —0.01285(38) (the experimental estimate of Ref. [jl|), r] = 0.0381(3), 
re = 1.96(2), rg = 1.40(15) and rio = ~13(7). This discrepancy is mainly due to the dif- 
ferent value of a, since the ratio A'^/A~ is particularly sensitive to it. This fact is also 
suggested by the phenomenological relation [p3] y4+/A^ 1 — ia. 
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TABLE XVII. Estimates of jA obtained in different approaches. 



IHT PR 


(i=3 exp. 


e-exp. 


experiments 


1.062(4) 


1.056(4) II 


1.029(13) H 


1.0442 igi 


1.055(3) |2|] 






1.067(3) g 








1.058(4) m 








1.088(7) |6| 



We observe a discrepancy with the experimental result reported in Refs. [|l],3, j A' = 
1.0442. However, we note that in the analysis of the experimental data of Ref. |^ the 
estimate of A'^/A'^ was strongly correlated to that of a; indeed A'^/A'^ was obtained by 
analyzing the high- and low-temperature data with a fixed to the value obtained from the 
low-temperature data alone. Therefore the slight discrepancy for A'^ /A' that we observe is 
again a direct consequence of the differences in the estimates of a. 

For the other universal amplitude ratios we quote as our final estimates the results 
obtained by using scheme (A) with n = 1: 



R+ = {A+y/''f+ = 0.355(3), 


(85) 


aA+C+ 
Rc = — ^ = 0.127(6), 


(86) 


^x= ^sc'^y =1-35(7), 


(87) 


^4 = = koP = 7.5(2), 


(88) 


= lim z-^F{z) = 0.0302(3), 

2— >00 


(89) 


c/ = 4(2). 


(90) 



These results are substantially equivalent to those reported in Ref. [26| 



V. THE TWO-POINT FUNCTION OF THE ORDER PARAMETER IN THE 

HIGH-TEMPERATURE PHASE 

The critical behavior of the two-point correlation function G{x) of the order parameter 
is relevant to the description of scattering phenomena with light and neutron sources. 

In the HT critical region, the two-point function G{x) shows a universal scaling behavior. 
For A;,m— ^0(m = l/,^ and ^ is the second-moment correlation length) with y = k'^/m'^ 
fixed, we can write fH^ 

g{y) = x/G{k). (91) 
The function g{y) has a regular expansion in powers of y: 

oo 

g{y) = l + y + Y,c^y\ (92) 

i=2 
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Two other quantities characterize the low- momentum behavior of g{y): they are given by 
the critical limit of the ratios 



SM = ml^p/m 



Sz = xm^/Z^ 



:ap5 



(93) 
(94) 



where mgap (the mass gap of the theory) and Zgap determine the long-distance behavior of 
the two-point function: 



47r|x| 



(95) 



If Uo is the negative zero of g{y) that is closest to the origin, then, in the critical limit, 
Sm = -yo and Sz = dg{y)/dy\y^y^. 

The coefficients Cj can be related to the critical limit of appropriate dimensionless ratios 
of spherical moments of G{x) [^5|,|l6| and can be computed by analyzing the corresponding 
HT series in the 0^ and in the dd-XY models, which we have calculated to 20th order. We 
report only our final estimates of C2 and C3, 



C2 = -3.99(4) X 10"^ 
C3 = 0.09(1) X 10-^ 



and the bound 



|C4| < 10"^ 

As already observed in Ref. ESI, the coefficients show the pattern 



\Ci\ < |Ci_i| < ... < IC2I < 1 



for 



i > 3. 



(96) 
(97) 



(98) 



(99) 



Therefore, a few terms of the expansion oi g{y) in powers of y provide a good approximation 
in a relatively large region around y = 0, larger than \y\ < 1. This is in agreement with the 
theoretical expectation that the singularity of g{y) closest to the origin is the three-particle 
cut (see, e.g., Refs. p8|J5P|, 
expansion of g{y) is Vg 
should have 



If this is the case, the convergence radius Vg of the Taylor 
9Sm- Since, as we shall see, Sm ~ 1, at least asymptotically we 



Ci+l 



1 

Cj. 

9 



(100) 



This behavior can be checked explicitly in the large- limit of the A^-vector model 
Assuming the pattern ([ 



we may estimate Sm and Sz from C2, C3, and C4. We obtain 



SaI = 1 + C2 - C3 + C4 + 2C2 + ... 

5"^ = 1 - 2c2 + 3c3 - 4c4 



2cl + ... 



(101) 
(102) 



where the ellipses indicate contributions that are negligible with respect to C4. In Ref. |^ 
the relation (|101|) has been confirmed by a direct analysis of the HT series of Sm- From 
Eqs. ([TUT| ) and we obtain 
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Sm = 0.999592(6), Sz = 1.000825(15). (103) 

These results improve those obtained in Ref. by using HT methods in the standard 



XY model and field-theoretic methods, such as the e expansion and the fixed-dimension g 
expansion. 

For large values of y, the function g{y) follows the Fisher-Langer law 
The coefficients have been computed in the e expansion to three loops |6^, obtaining 



Ai ^ 0.92, A2 ^ 1.8, A3 ^ -2.7. (105) 

In order to obtain an interpolation that is valid for all values of y, we will use a phenomeno- 
logical function proposed by Bray |6^. This approximation has the exact large-y behavior 



and its expansion for ?/ ^ satisfies Eq. (|100|) . It requires the values of the exponents u 



a, and rj, and the sum of the coefficients A2 + A3. For the exponents we use of course the 
estimates obtained in this paper, while the coefficient A2 + A3 is fixed using the e-expansion 
prediction A2 + A3 = —0.9. Bray's phenomenological function predicts then the constants 
Ai and q. We obtain: 

Ai ^ 0.915, A2 ^ -24.7, A3 ^ 23.8, (106) 

C2 ^ -4.4 X 10-^ C3 ^ 1.1 X 10-^ C4 ^ -5 X 10"^. (107) 

The results for Ai, C2, C3, and C4 are in good agreement with the above- reported estimates, 
while A2 and ^3 differ significantly from the e-expansion results ( |105[ ). Notice, however, 
that, since |a| is very small, the relevant quantity in Eq. (|104|) is the sum A2 + A3 which is. 



by construction, equal in Bray's approximation and in the e expansion: in other words, the 
function (|104|) does not change significantly if we use Eq. (|105|) or Eq. (|106|) for A2 and A3. 



Thus, Bray's approximation provides a good interpolation both in the large-y and small-?/ 
regions. 

APPENDIX A: THE MONTE CARLO SIMULATION 

1. The Monte Carlo algorithm 

At present the best algorithm to simulate iV-vector systems is the cluster algorithm 
proposed by Wolff |71| (see Ref. |^ for a general discussion). However, the cluster update 



changes only the angle of the field. Therefore, following Brower and Tamayo [[7^, we add a 
local update that changes also the modulus of the field. 

a. The wall- cluster update 



We use the embedding algorithm proposed by Wolff |71] with two major differences. 



First, we do not choose an arbitrary direction, but we consider changes of the signs of the 
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first and of tlie second component of the fields separately. Second, we do not use the single- 
cluster algorithm to update the embedded model, but the wall-cluster variant proposed in 
Ref. [^. In the wall-cluster update, one flips at the same time all clusters that intersect a 
plane of the lattice. In Ref. we found for the 3D Ising model a small gain in performance 
compared with the single-cluster algorithm. 

Note that, since the cluster update does not change the modulus of the field, identical 
routines can be used for the 6'^ model and for the dd-XY model. 



b. The local update of the (j)'^ model 

We sweep through the lattice with a local updating scheme. At each site we perform a 
Metropolis step, followed by an overrelaxation step and by a second Metropolis step. 
In the Metropolis update, a proposal for a new field at site x is generated by 

0f) = 0(1)+ c(r«- 0.5), 

0f = 0i'^ + c(r(^)-O.5), (Al) 

where r^^) and r^^) are random numbers that are uniformly distributed in [0, 1). The proposal 
is accepted with probability 

A = min[l, exp{-H' + H)]. (A2) 

The step size c is adjusted so that the acceptance rate is approximately 1/2. 
The overrelaxation step is given by 




— * — * 

where 0„ = "^y^nnix) nn(a;) is the set of the nearest neighbors of x. Note that this 

step takes very little CPU time. Therefore, it is likely that its benefit out-balances the CPU 
cost. 



c. The local update of the dd-XY model 

We sweep through the lattice with a local updating scheme, performing at each site one 
Metropolis update followed by the overrelaxation update (|A3|). 

In the Metropolis update, the proposal for the field (px at the site x is given by 

0L = (O,O) for 10,1 = 1 

$'x — (cos(q;), sin(a)) for |0a,| = 0, (A4) 

where a is a random number with a uniform distribution in [0, 2tt). This proposal is accepted 
with probability 

A{$'^, 4) = min[l, expi-H' + H)] = min[l, exp(/5 4 " (0^ - 4) + ^(0^' - $1))], (A5) 

where 0„ = J2yenn{x) nearest-neighbor spins. We will prove in Sec. |A 1 e| 

that this update leaves the Boltzmann distribution invariant. 
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d. The update cycle 

Finally, we summarize the complete update cycle: 
local update sweep; 

global field rotation, in which the angle is taken from a uniform distribution in [0, 27r); 
6 wall-cluster updates. 

The sequence of the 6 wall-cluster updates is given by the wall in 1-2, 1-3 and 2-3 plane. 
In each of the three cases, we update separately each component of the field. 



e. Proof of the stationarity of the local update of the dd-XY model 

Stationarity means that the update leaves the Boltzmann distribution invariant. Since 
our update is local, it is sufficient to consider the conditional distribution of a single spin 
for given neighbors 

PbM = \ exp(/? 4>x-$n + D^l), (A6) 

where 

z = j dn{(j)^) exp(/? (j)x-(l>n + D$l) (A7) 
and dfi{(j)ri.) is defined in Eq. (^. We have to prove that 

n(C) = / w{$'^, 4) (A8) 



is satisfied, where W is the update probability defined by Eqs. ( |A^ ) and ([A5| ) . 
Using Eq. (H), the right-hand side of Eq. (|A8|) can be rewritten as 



W{$',, 0) Pb(0) + J dn{(j)^) W{$',, Pb(0,), (A9) 



where dQ{(f)x) is the normalized measure on the unit circle. 
Let us first consider the case \d>'\ = 1. Then, we have 



and 



iy(0l.,O)n(O) = ^A(0^,O), (AlO) 



= n(0;)-^A(0;,o), (All) 



where we have used the property 
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A{0,ij)P,{i:)=A{ij,0)P,{0), (A12) 



valid for \ip\ = 1. Summing the two terms we obtain Eq. ( |A8| ) as required. 
For {(h'l = 0, we have 



W{0, O)Pfe(O) = i (^1 - y dn{ip) A{^p, 0)^ , (A13) 



dn{<f),) w{o, 0,) pS^) = - I dn{tij) A{tij, 0) (am) 



and 



where we have used again Eq. ( [A12| ). Summing the two terms we obtain Eq. ([A8|) as required. 



2. Measuring Za/Zp 



One of the phenomenological couphngs that we have studied is the ratio Za/Zp of the 
partition function Za of a system with anti-periodic boundary conditions (a.b.c.) in one 
direction and the partition Zp with periodic boundary conditions (p.b.c.) in all three direc- 
tions. A.b.c. are obtained by multiplying the term (I)x4'y in the Hamiltonian by —1 for all 
X = (Li, X2, xs) and y = (1, X2, x^). This ratio can be obtained using the so-called boundary- 
flip algorithm, applied in Ref. |Q to the Ising model and generalized in Ref. |44| to general 
0(iV)-invariant non-linear a models. 

In the boundary-flip algorithm, one considers fluctuating boundary conditions, i.e., a 
model with partition function 



^fluct — Zp + Za — 

Jb=±l 



D[0] exp 



P'^J{xy)4>x(py + 



(A15) 



where Ji^^y) = Jb for ^ = {Li, X2, x^) and y = (1, X2, x^), and Ji^^y) 
Jh = —1 correspond to p.b.c. and a.b.c. respectively. 

In this notation, the ratio of partition functions is given by 



1 otherwise. Jb = ^ and 



(A16) 



where the expectation value is taken with fluctuating boundary conditions. 

In order to simulate these boundary conditions, we need an algorithm that easily allows 
flips of Jb- This can be done with a special version of the cluster algorithm. For both 
components of the field we perform the freeze (delete) operation for the links with probability 



p, = min[l,exp(-2/?J^.,)0(f)0j'))], 



(A17) 



where ^i^"* is the chosen component of (f)x- The sign of Jb can be flipped if there exists, for the 
first as well as the second component of the field, no loop of frozen links with odd winding 
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number in the first direction. In Ref. it is discussed how this can be implemented. For 
a more formal and general discussion, see Ref. 



75 . Note that for Jh 



—1 the flip can 

always be performed. Hence, as F. Gliozzi and A. Sokal have remarked [^], the boundary 
flip needs not to be performed in order to determine Za/Zp. It is sufficient to use p.b.c. and 
check if the flip to a.b.c. is possible. Setting 6 = 1 if the boundary can be flipped and 6 = 
otherwise, we have 



Za 



(A18) 



where the expectation value is taken with p.b.c. 



3. Checks of the program 



a. Schwinger-Dyson equations 



The properties of the integration measure allow to derive an infinite set of nontrivial 
equations among observables of the model. Here we have used two such equations to test 
the correctness of the programs and the reliability of the random-number generator. For a 
more general discussion of such tests, see Ref. [[T^ . 

For the 0'^ model, the partition function remains unchanged when, at site x, the first 
component of the field is shifted by ip: (j)^^ — > ^i^"* + ip. We obtain, using also the 0(2) 
invariance 



1 + ij)^ 



Z 



dtp 



2 



+ 1, 



(A19) 



where {■■■)^ indicates that the Boltzmann factor is taken with the shifted field at the site x. 
A second equation, which is valid for the (p^ model as well as for the dd-XY model, can be 
derived from the invariance of the measure under rotations: 



0f 



cosl^a^ 
sinfa) 



(1) _ 



sm a 



^(2) 



0« + cos(a) 0f 



Taking the second derivative of the partition function with respect to a yields 



1 d^Z(a) 



Z da^ 



a=0 




y£nn(x) 



■'xfyli 



(A20) 



(A21) 



where y e nn(a;) indicates that the sum runs over the six nearest neighbors of x. 

We checked Eq. ( [A19| ) for all our simulations of the 0^ model. In the simulation we 
averaged Eq. (|A19|) over all sites of the lattice to reduce the error bars. For most of the sim- 
ulations the right-hand side of Eq. ( [A19| ) differed by zero by less than 2 standard deviations. 
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Only in one case (L = 8, A = 2.2) the difference was about 3 standard deviations. The 



weighted average of the right-hand side of Eq. ( A19 ) over all our simulations is —8(9) x 10'^. 
Hence, there is no indication for a program error or a problem with the random-number 
generator. 



We implemented Eq. (|A21| ) in all simulations of the standard XY model and, unfortu- 



nately, only in the most recent simulations of the dd-XY model. For L = 96, D = 1.02, and 
P = 0.56379 we found Pm = 0.5637896(24) from 575,000 measurements, where 
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We have measured jSm in all our simulations of the XY model. In most cases the deviation 
of (3rn from f3 = 0.454165 was less than one standard deviation. The largest deviation 
was — (3 = —0.0000049(22) for L = 28. The weighted average over all simulations is 
(3^ -(3 = -0.00000024(52). 

Also this check does not indicate a program error or a problem with the random-number 
generator. 

b. Checks of the Taylor expansion 

As a test of the MC program and of the analysis software, we simulated the 0^ model 
for L = 4 and A = 2.1 at the following values oi 13: /3 = 0.485, 0.490, 0.495, 0.500, 0.505, 
0.510, 0.515, and 0^520. 

We computed R with Rij = {Za/Zp)f = 0.3202 and Rij = (6nd/^)/ = 0.5925 and 
i?2 = and i?2 = Uq for all these simulations, using a third-order Taylor expansion. The 
results for R are summarized in Table PCVII1| . They show that there is a large interval in 



which the method works: indeed, the results for R for (3s = 0.505, 0.510 and 0.515 agree 
within two standard deviations, although the variation of U4 and Uq at Ps is several hundred 
standard deviations. In addition we have gained information about the range of (3 where 
the extrapolation works with the desired accuracy: 

1/3, -(3f\< 0.005 X (L/4)i/^ (A23) 

The factor {1/4^^" takes care of the fact that the slope of the couplings R scales like L^^'^. We 
have carefully checked that this requirement is always fulfilled in our simulations. Therefore, 
we are confident that the extrapolation in (3, using the Taylor expansion, is implemented 
correctly. 



4. CPU-time, acceptance rate, cluster size, and autocorrelation times 

We have used ANSI C to implement our simulation programs. We have used our own 
implementation of the G05CAF random-number generator from the NAG-library. The 
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TABLE XVIII. Test of the Taylor expansion. Simulations with L = A, (5f ^ 0.50773 for 
{Za/Zp)j, and j3f ~ 0.50994 for (^2nd/^)/- In the first column we give the value [3s where the 
simulations have been performed (i.e., the configurations have been generated with a weight pro- 
portional to the Boltzmann factor that corresponds to (ig)- In the columns 2 up to 5 we give the 
results for four choices of R. These choices are labelled by R2 at Rij- Finally in columns 6 and 7 
we give C/4 and Uq at f3s for comparison. 



Ps Ua at {Za/Zp)f Uq at {Za/Zp)f Uj at (6nd/-^)/ U& at (6nd/^)/ Ua at 13s 



Uq at f3s 



0.485 
0.490 
0.495 
0.500 
0.505 
0.510 
0.515 
0.520 



1.249912(46) 
1.249589(42) 
1.249431(40) 
1.249343(39) 
1.249365(39) 
1.249458(40) 
1.249373(42) 
1.249313(45) 



1.76951(17) 
1.76749(15) 
1.76671(14) 
1.76634(14) 
1.76642(13) 
1.76678(14) 
1.76641(14) 
1.76616(15) 



1.239726(58) 
1.239863(51) 
1.239866(46) 
1.239805(43) 
1.239824(41) 
1.239885(40) 
1.239802(39) 
1.239804(39) 



1.73513(21) 
1.73434(18) 
1.73396(16) 
1.73364(15) 
1.73370(14) 
1.73394(13) 
1.73359(13) 
1.73358(13) 



1.360445(71) 
1.334375(67) 
1.309184(63) 
1.284661(60) 
1.261574(57) 
1.239660(54) 
1.218859(50) 
1.199535(47) 



2.17525(28) 
2.07444(25) 
1.98000(23) 
1.89063(22) 
1.80882(20) 
1.73318(18) 
1.66302(17) 
1.59940(15) 



G05CAF is a linear congruential random-number generator with modulus m = 2^^, mul- 
tiplier a = 13^^ and increment c = 0. Most of our simulations have been performed on 450 
MHz Pentium III PCs running the Linux operating system. 

For D — 1.03, the average size of the clusters per volume decreases from about 0.52 for 
L = 5 to about 0.14 for L = 80. For large lattices the behavior of the size of the wall clusters 
per volume is roughly given by 1.198 x /.-O-^ss, Por A = 2.1 we get 1.276 x /.-O-^ss^ 

The acceptance rate of the local update of the dd-XY model at D = 1.03 is about 0.273 
for L — 5. It increases slowly and reaches about 0.309 for L — 80. At the same time the 
density oi cf) — (0, 0) spins increases from about 0.150 to 0.16628. 

The acceptance rate of the local update of the (j)^ model at A = 2.1 is about 0.395 for 
L = 5 and increases to 0.4098 for L = 80. We have chosen a step size c = 2.0 throughout. 

For the dd-XY model for 16 < L < 80 one cycle of the update plus the measurement 
takes about x 3.3 x 10~^ seconds on a 450 MHz Pentium HI CPU. In this range of lattice 
sizes the decrease of the wall-cluster size per volume and the increase of memory access 
times seems to cancel almost exactly (by chance). 

Roughly 2/3 of the CPU-time is spent with the cluster update (which incorporates the 
measurements needed for Za/Zp). Again about 2/3 of the remaining time is spent in the 
local update. 

The CPU time required for the update of the cf)^ model is shghtly larger due to a different 
local update and due to slightly larger wall clusters than in the dd-XY model. 

In total about 5 years on a single 450 MHz Pentium HI CPU were used for our study. 
(We had 8 CPUs available for the study.) 

We performed additional runs of the 0^ model at A = 2.1 and P — 0.50915 and the dd- 
XY model at = 1.03 and /3 = 0.5628, where we had stored all measurements of the square 
of the local field, the energy, the magnetization, F, and the boundary variable. Hence, we 
could compute autocorrelation functions and times from these data. We simulated lattices of 
size L = 8 up to L = 64 for this purpose. In all cases we performed 200,000 measurements. 
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TABLE XIX. Integrated autocorrelation times for the (j)'^ model at A = 2.1 and /? = 0.50915, 
and the dd-XY model at D = 1.03 and P = 0.5628. We have always performed 200,000 cycles. 
The autocorrelation times are given in units of update cycles, te is the integrated autocorrelation 
time of "^(^xy) 4>x4>y, Tcffi of Ylix '^x^ '^x magnetic susceptibility, Tp of the Fourier transform of 

the field at minimal nonvanishing momentum, and Tb of the boundary variable, "trunc" gives the 
time at which the summation of the autocorrelation function was truncated. Note that trunc was 
chosen to be the same for all observables. 



L 


TE 






TF 


n 


trunc 


(j)'^ model 


8 


2.61(5) 


2.36(4) 


2.30(4) 


0.79(1) 


1.67(3) 


15 


12 


3.12(6) 


2.72(5) 


2.68(5) 


0.86(2) 


1.90(4) 


18 


16 


3.39(7) 


2.85(6) 


2.81(6) 


0.89(2) 


1.99(4) 


20 


24 


4.06(9) 


3.25(7) 


3.31(7) 


0.94(2) 


2.29(5) 


25 


32 


4.54(11) 


3.51(9) 


3.57(9) 


0.94(2) 


2.53(6) 


30 


48 


5.65(15) 


4.22(11) 


4.23(11) 


1.01(2) 


2.93(8) 


35 


64 


5.89(17) 


4.36(12) 


4.26(12) 


1.01(3) 


3.08(9) 


40 


dd-XY model 


8 


2.80(4) 


1.78(3) 


2.47(4) 


0.78(1) 


1.79(3) 


15 


12 


3.21(6) 


2.01(4) 


2.74(5) 


0.81(2) 


1.91(3) 


18 


16 


3.57(7) 


2.25(5) 


2.97(5) 


0.88(2) 


2.13(4) 


20 


24 


4.23(10) 


2.66(6) 


3.44(8) 


0.93(2) 


2.42(5) 


25 


32 


4.92(12) 


3.14(7) 


3.83(10) 


0.99(2) 


2.75(6) 


30 


48 


5.81(15) 


3.65(10) 


4.17(11) 


0.97(3) 


2.95(8) 


35 


64 


6.58(19) 


4.24(12) 


4.55(13) 


1.01(3) 


3.22(9) 
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A measurement was performed after each update cycle. 

Our results for the integrated autocorrelation times are summarized in Table pClX| . All 
autocorrelation times are given in units of measurements. We see that the autocorrelation 
times slightly increase with increasing L. At a given lattice size, the autocorrelation times 
for the (f)^ model are a little smaller than those for the dd-XY model. Likely, this is caused 
by the particular implementation of the local updates. Among the observables that we have 
studied, the energy has the largest autocorrelation time. It is interesting to note that the 
autocorrelation time of (l)x shows no sign of severe critical slowing down, despite the 

^ 2 

fact that (px is only changed by the local update. 



APPENDIX B: ANALYSIS OF THE HIGH-TEMPERATURE EXPANSIONS 

In this appendix we report a detailed discussion of our HT analyses. This detailed 
description should allow the reader to understand how we determined our estimates and the 
reliability of the errors we report, which are to some extent subjective. 
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1. Definitions and HT series 



Using the linked-cluster expansion technique, we computed the 20th-order HT expansion 
of the magnetic susceptibihty 

X = Y.^<Pama{x)), (Bl) 

X 

of the second moment of the two-point function 

m2 = ^x2(0„(O)0,(x)), (B2) 

X 

and therefore, of the second-moment correlation length 

5^ = p. (B3) 

Moreover, we calculated the HT expansion of the zero-momentum connected 2j-point 
Green's functions X2j 

X2j= (<^"i(O)<^«i(^2)--0a,(2;2i-l)0a,(a;2i))c (B4) 

X2,...,X2j 

ix = X2)- More precisely, we computed X4 to 18th order, xe to 17th order, xs to 16th order, 
and Xio to 15th order. The series for the 0^ Hamiltonian with A = 2.07 and the dd-XY model 
with D = 1.02 are reported in Tables pCX| and pCXl| . The HT series of the zero-momentum 
four-point coupling (74 and of the coefficients that parametrize the equation of state can 
be computed using their definitions in terms of X2j and i.e.. 



3iV X4 .R,. 



and 



in 5(iV + 2) x6X2 



3 ' 



= 980 - 280(iV + 2) X6X2 35{N + 2f xsxl 

3(iV + 4) xl 9{N + 4){N + 6) xl 
= 15400 - 7700(iV + 2) XeX2 350(iV + 2)' xlxl 



+ 



(iV + 4) xl (iV + 4)2 x\ 
1400(iV + 2)2 xsxl 35(iV + 2)3 xioxl 



3(Ar + 4)(Ar + 6) xl 3(Ar + 4)(Ar + 6)(Ar + 8) xl 
The formulae relevant for the XY universality class are obtained setting N = 2. 
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TABLE XX. CoefHcients of the HT expansion of m2, x, X4, Xe^) XSj and xio- They have been 
computed using the cf)^ Hamiltonian with A = 2.07. 



z 


M2 


X2 


X4 








0.82195468340525626553069 


-0.18234682673209145113698 


1 


1.01341425235775258955047 


2.02682850471550517910095 


-1.79856993883835218414950 


2 


4.99788354573054604609481 


4.39836023278442865137831 


- 10. 1455394390643744635266 


3 


17.0521428795050699809477 


9.45608303384639636168661 


-44.4379274874486579169902 


4 


50.2802232377961742954821 


19.7509483421835061143003 


-166.911247402827967481247 


5 


136.081955345771888598704 


41.0953224517675210993058 


-566.231817578861031037728 


6 


349.014986520228913452360 


84.3784681474012731739965 


- 1783.81074566131415934408 


7 


861.072204516234675501406 


172.831419974338990927678 


-5314.73425458572561425498 


8 


2065.11115559310738635122 


351.412467687273478895617 


-15153.8054671243945401658 


9 


4843.65801296852863958594 


713.327141375620774069366 


-41706.6393990107185189143 


10 


11163.1410843891753269547 


1441.29537543678062498269 


-111480.350225534613432415 


11 


25357.2828059634753398090 


2908.62198285250042113961 


-290779.553330813045449802 


12 


56914.8160479537305800002 


5851.22700855642135604880 


-742792.316757751368697227 


13 


126448.957588634753901037 


11759.8313973078166535503 


— 1863674.70918289087066236 


1 A 


07sc;n/i 7Q/iQQS07nofinc;i 1 o/i/i 

Z 1 00U4I:. < y4oooZ ( UZDU0iiZ44 


zooou.DyoyzoooaZ4iUDyooDo 


— 4DUooio.yizzyoio4oo < Z4iy < 


15 


608 ^ /5.0440o8biy8o4y011z4 


4 r z4o.yziiyUzoiooboz / /4io 


— iizi4y4o.zoioo4oozzyoo5yb 


lb 


13ziy48.z7b88188y4433y5o3 


y4o08. 3475046413360663308 


—2byy 1439. 1725436301683188 


17 


28538z3.y4y33b432z0583008 


188y24.8113y7083165bb0ybl 


—64258568.1525362515312215 


1 o 
lo 


oizoyoU.ozUUooobozi rozoz4 


377150. olDZzo4yz7oyyo3104 


— iOi4y2 roU.iU42iOOUUDb4Uo2 


ly 


ioiUi4D / .OODzyozyzUoD / ozi 


( ozoo4. 1 zooDo4oUozi iyy4y i 




on 
20 


z7ooyizy.D7DlDz7014bo7zb4 


i4yyoyo. io0i4 / oUdzoU4o4Uz 




i 


X6 


X8 


Xio 





0.41197816889670188853628 


-2.04244608438484921752518 


17.4966181373390399118217 


1 


8.09031485009238719234751 


-65.5324533503954629382215 


827.907846284122870554505 


2 


80.5579033480687931078788 


-993.099544110386544944672 


17583.1558423294442666385 


3 


569.040955700888943670354 


-10125.9430216092445934264 


242162.737636025061958065 


4 


3235.75380473965447264355 


-79870.4410385967305759821 


2508010.66864972545742050 


5 


15823.8724636340283846359 


-524859.679867495864300189 


21148750.0733186453907577 


6 


69189.8243270623365490494 


-3005237.58408599956662949 


152451568.855786252573448 


7 


277430.904691458150896928 


-15443601.8321215251366505 


970471959.824139222337835 


8 


1038008.04844006612521747 


-72717890.9115560213984880 


5582186363.05878798571648 


9 


3669720.02487068188513696 


-318510657.462768013289640 


29507970388.4631661994252 


1 n 


IZO/^ioZ. 1 DUZUDoOOUUiyOUo 


1 "^1 9(^fiR7f^7 1 70/1 ^^^^7^1 *^QRQ 

— ioiZDOo ( D ( . 1 1 u^ooo 1 Dioyoy 


1 /L^9n'^/1RQ99Q zlfil '^0000^077 


11 


40084338.4384102309965122 


-5135470803.97668705323439 


671866031736.491009927361 


12 


125446404.887628912440242 


- 19206423870.3058981804587 


2946715041148.83369711661 


13 


381003987. 146313729964559 


-69057827061.9724837730205 


12330196249913.8306690405 


14 


1127156412.88952101681812 


-239824993302.080594778804 


49488831955055.0830738902 


15 


3257906807.50792443987963 


-807540440278.527392223891 


191378645936343.575972510 


16 


9223412391.97191941389912 


-2645013087720.90303565763 




17 


25631282620.7774958190658 
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TABLE XXI. Coefficients of tfie HT expansion of m2, Xi X4) X6)5 XSj and xio- They have been 
computed using the dd-XY Hamiltonian with D = 1.02. 



I 


M2 


X2 


X4 








0.73497259946651881066155 


-0.12952381667498436104661 


1 


0.81027708294985783008938 


1.62055416589971566017876 


- 1 . 14235747481325709236407 


2 


3.57318872366283058263222 


3. 19240289872507821851086 


-5.86566452750871250322137 


3 


1 1 .0006482367650529136739 


6.24412164584056694746990 


-23.4964691362678107734860 


4 


29.3868020693268986850123 


11.8078991080110710061038 


-80.6562687454557047131467 


5 


72.0537292098720051000559 


22.2452167397051997713027 


-249.761461812242731292855 


6 


167.391323293772704647731 


41.3083108337880277068933 


-717.220908333053642342640 


7 


373.956365385479749391430 


76.5360051001560916440577 


-1945.68221342894368675090 


8 


811.915341705911027278243 


140.6796608851 19763690958 


-5046.15791369304901147071 


9 


1723.52034635335232533650 


258.178251744079668880881 


- 12622.3575801936426660898 


10 


3594.34333565606104441219 


471 .478755943757026593209 


-30642.5601039304090869577 


11 


7386.64345245316356464020 


860.010504587591344158158 


- 72548.9796473949259402584 


12 


14997.4679478666032080397 


1563.48753262138558216805 


-168135.618026803729643010 


13 


30136.9726456923668836816 


2839.86987807993137143111 


—382565.790192788611922604 


1 /I 


Duuzy . ioo ( oyoyozzu4Ui iooo 


0i40.o40UoDooo ( 4 < UoDzyiy4 


— oODDzy.y4oOoUOyOD4DZ4Dyoz 


15 


iioDOD.OiDy5DZD/OZ / iOOZZO 


931 ^^.6 1 z06o91z4083z59d1zd 


— ioyio5i.24U ( 6y 1 2ioU4yU224 




z6Zy /9. 69980 ( 4540ol09o5zy 


16841.06600 r 613018 1 494565 


— 4124106. iiioU66ooyy4i2oo6 


17 


454746.664171304150538747 


304zl.55731673384658058z5 


— 88y3532.37374656560556y00 


lo 


oooft^?n ftOyiTft/i cco/( /( 1 noi OftCCO 
oozyoU.yz4 r y4oo44iUoizyoo 


CC/IQ7CC /( 70ftOftn^i 1 o 1 n^i ccono^ift 
O4o7o.47zyo90613106o30o69 


1 Oft07ftCCO ft^?ftni CCyl /( OftCCftI ooo 

— Ioyo7yo3.y6yoio44oyoyi3o3 


1 ft 


1 Tn^iQQn A'7nft'70'7^^/( c^qcqqi nft 
i / UoooU.O /UU ( z / 04oooooiUU 


nOOQC nQ7ni f?Cft7nQ^^C^Q71 QQQ 

yoydo.yo < Uiooy < Uoooo t looo 




on 
20 


3zooo69.770zo650z4zo4oz76 


1 701 oo nncc7ccn^ini ftnc^cc7nft7^? 
178182.095700601900070976 




i 


X6 


Xs 


Xio 





0.19923804166181643967757 


-0.67796674603175173967960 


4.06523419917732615691360 


1 


3.64240617023376584468207 


-20.6297276122379010898516 


182.870340299506381670926 


2 


33.7843552454418076210054 


-295.585300853935935335485 


3697.42581896510915833415 


3 


221.913477156277495814452 


-2838.14010017098792994188 


48330.0439419608161526720 


4 


1169.53606787476515583468 


-20984.8625853355842977237 


472993.887593911146392516 


5 


5284.11768819259530573126 


-128723.340916169008103149 


3752687. 19462636307477999 


6 


21286.1074543663282451107 


-685456.078166050751079896 


25350793.0882363704732840 


7 


78446.3667728434480782029 


-3265581.01409539304748618 


150699914.609020169056435 


8 


269224.069430361947223693 


-14216412.9210915918573320 


806981988.823436545704567 


9 


871590.289814654747926984 


-57438991.7197952962871671 


3960575199.72383885451240 


1 n 


ZDo / 4 / O.OUoy4U 1 o4il lOUUOO 


91 7Q971 1 1 RR'^RR7'^*^fin*^n'^n 


loUOZoloOZo.UO / DO'^z^^OU^OO 


11 


7951104.44962467048194209 


-783550202.710253766125616 


77211058033.5079429953492 


12 


22703108.4838172682261707 


-2689172405.44915075129067 


312454827813.125173762158 


13 


62855482.2396162632890035 


-8861507312.08886862026490 


1204410146707.71423350686 


14 


169374342.873715022154662 


-28171979758.2629139055303 


4446806950469.81858484485 


15 


445613055.420361814905254 


-86751935057.8910676591399 


15798640107921.6847945606 


16 


1147647339.61141521962212 


-259625721060.742527022681 




17 


2899724764.50550187290555 
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TABLE XXII. MC estimates of A 



/3c 


A = 2.00 


0.5099049(15) 


A = 2.07 


0.5093853(24) 


A = 2.10 


0.5091507(13) 


A = 2.20 


0.5083366(24) 


D = 0.90 


0.5764582(24) 


D = 1.02 


0.5637972(21) 


D = 1.03 


0.5627975(14) 


D = 1.20 


0.5470377(26) 



2. Critical exponents 

In order to determine the critical exponents 7 and v from the HT series of x ^^nd j3 
respectively, we used quasi-diagonal first-, second-, and third-order integral approximants 
(lAl's, IA2's and IA3's respectively). Since the most precise results are obtained by using 
the MC estimates of /3c to bias the approximants, we shall report only the results of the 



biased analyses. The values of f5c used in the analyses are reported in Table PCXII . 

Given an nth-order series /(/3) = 'Y^=o^iP^i its fcth-order integral approximant 
[mfc/mfc_i/.../mo//] lAA; is a solution of the inhomogeneous /cth-order linear differential 
equation 

pmf^'\fi) + Pk-mf'~-'\p) + ... + pmf^'\P) + pmiiP) + r{p) = o, (b7) 

where the functions Pi{f3) and R{l3) are polynomials of order rrii and I respectively, which 
are determined by the known nth-order small-/? expansion of /(/3) (see, e.g., Ref. [p!0|]). 
We consider three types of biased lA/c's: 

(i) The first kind of biased IAA;'s, which will be denoted by blA/c's, is obtained by setting 

Pfc(/3) = (l-/3//5,)p,(/3), (B8) 

where Pk{(3) is a polynomial of order — 1. 

(ii) Since on bipartite lattices (3 = —(3c is also a singular point associated to the antifer- 
romagnetic critical behavior [^, we consider lA/c's with 



PM = {l-^'/Pl)pM, (B9) 

where Pk{(3) is a polynomial of order — 2. We will denote them by h±lAk's. 

(iii) Following Fisher and Chen [|l9l, we also consider lAfc's where the polynomial asso- 
ciated with the highest derivative of f{f3) is even, i.e., it is a polynomial in (3"^. In this case 
ruk is the order of the polynomial Pk as a function of /5^, i.e., Pk = X]j=o ^i'^^"'- Thus, in 
order to bias the singularity at Pc, we write 

Pk{f3) = {l-P'/P!)pk{P'), (BIO) 
where Pk{(3) is a polynomial in /3^ of order rUk — 1. We will denote them by bFCIA/c's. 
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In our analyses we consider diagonal or quasi-diagonal approximants, since they are 
expected to give the most accurate results. Below, we give the rules we used to select the 
quasi- diagonal approximants. We introduce a parameter q that determines the degree of 
off-diagonahty allowed (see below) . In order to check the stability of the results with respect 
to the order of the series, we also perform analyses in which we average over the results 
obtained with series of different length. For this purpose, we introduce a parameter p and 
perform analyses in which we use all approximants obtained from series of n terms with 
n > n > n — p. 

We consider the following sets of lA/c's: 

(a) [mi/mo/k] blAl's with 

n > mi + rriQ + k + l>n — p, 

Max [[{n - 1)/3J -q,3]< mi, mo, k<\{n- l)/3] + q. (Bll) 

(b) [mi/mo/A;] b±IAl's with 

n > nil + nio + k > n — p, 

Max [[{n - 1)/3J - g, 3] < mi, mo, k<\{n- l)/3l + q. (B12) 

(c) [mi/mo/k] bFCIAl's with 

n > mi -|- mo + k + l>n — p, 

Max[[(n - 1)/3J - q,3] <mi,mQ,k< [(n-l)/3] + q. (B13) 

(d) [m2/mi/mo/k] bIA2's with 

n > 1712 + mi -|- mo -|-A;-|-3>n — p, 

Max[[(n-3)/4j - q,2] < m2 - l,mi,mo,k < [(n-3)/4] + q. (B14) 

(e) [m2/mi/mo/k] b±IA2's with 

n > 1712 + mi + mo + k + 2>n — p, 

Max[[(n-3)/4j - q,2] < m2 - 2,mi,mo, k < f(n-3)/4] + q. (B15) 

(f) [ma/mi/mo/A;] bFCIA2's with 

n > 1712 + mi -|- mo + k + 3>n — p, 

Max [ [(n - 3) /4J - 2] < ms - 1, mi, mo, A; < \{n - 3) /4] -f- q. (B16) 

(g) [1713/1712/1711/1710/ k] blAS's with 

n > 1713 + i7i2 -|- mi -|- mo -|-A;-|-5>n — p. 

Max [[(n - 5) /5J - g, 2] < mg - 1, 1712, mi, mo, k< \{n- 5) /5] + q. (B17) 

In the following we fix g = 3 for the lAl's and q = 2 for the IA2's and IA3's. 
For each set of lA/c's we calculate the average of the values corresponding to all nondefec- 
tive lA/c's listed above. Approximants are considered defective when they have singularities 
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close to the real (3 axis near the critical point. More precisely, we consider defective those 
approximants that have singularities in the rectangle 



Xmin < ^eP/Pc < a^max, |Im/5//3c| < 2/max- (B18) 

The values of Xmin, a^max and z/max are fixed essentially by stability criteria, and may differ 
in the various analyses. One should always check that the results depend very little on the 
chosen values of Xmin, x^ax, and z/max, by varying them within a reasonable and rather wide 
range of values. The domain (P18|) cannot be too large, otherwise only few approximants are 
left. In this case the analysis would be less robust and therefore less reliable. We introduce 
a parameter s such that 

•^min 1 'S, a^max 1 ~l~ S, 2/max ■S, (B19) 

and we present results for various values of s. We also discard some nondefective lA's — 
we call them outliers — whose results are far from the average of the other approximants. 
Such approximants are eliminated algorithmically: first, we compute the average A and the 
standard deviation a of the results using all nondefective lA's. Then, we discard those lA's 
whose results differ by more than n„a from A with = 2. We repeat the procedure on 
the remaining lA's, by calculating the new A and a, but now eliminating the lA's whose 
results differ by more than n^a with n„ = 3. The procedure is again repeated, increasing 
Ua by one at each step. This procedure converges rapidly and, as we shall see, the outliers 
so determined are always a very small part of the selected nondefective lA's. 

In the Tables pCXllI| and pCXlV| , we present the results for the critical exponents 7 and u 



respectively, obtained from the HT analysis of the 0^ Hamiltonian. In the Tables XXV and 



pCXVi] we report the results for the HT analysis of the dd-XY model. There we also quote the 
"approximant ratio" Va = {g — f)/ 1, where t is the total number of approximants in the given 
set, g is the number of nondefective approximants, and / is the number of outliers which are 
discarded using the above-presented algorithm; g — f is the number of "good" approximants 
used in the analysis; notice that g ^ f, and g — f is never too small. For each analysis, 
beside the corresponding estimate, we report two numbers. The number in parentheses, 
Ci, is basically the spread of the approximants for j3c fixed at the MC estimate. It is the 
standard deviation of the results obtained from all "good" lA's divided by the square root 
of Ta, i.e., ei = a/y/Va. Such a definition of ei is useful to compare results obtained from 
different subsets of approximants of the same type, obtained imposing different constraints. 
The number in brackets, 62, is related to the uncertainty on the value of Pc and it is estimated 
by varying in the range [Pc - ^Pc, Pc + ^Pd- 

The results of the analyses are quite stable: all sets of lA's give substantially consistent 
results. The comparison of the results obtained using all available terms of the series with 
those using less terms (in the Tables the number of terms is indicated explicitly when it 
is smaller than the number of available terms) and those obtained for p = 3 (i.e., using 
n, n — 1, n — 2, and n — 3 terms in the series) shows that the results are also stable with 
respect to the order of the HT series. Therefore, we do not need to perform problematic 
extrapolations in the number of terms, or rely on phenomenological arguments, typically 
based on other models, suggesting when the number of terms is sufficient to provide a 
reliable estimate. 
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TABLE XXIII. Results for 7 obtained from the analysis of the 20th-order HT series of x for 
the (j)^ Hamiltonian. The number n of terms used in the analysis is indicated explicitly when it is 



smaller than the number of available terms (n 


= 20). p = when its value 


is not explicitly given. 


A approximants 


ra 


7 


2.00 bIAl,=i/2 


(28 - 2)/48 


1.31755(11)[19] 


bIA2,=i/2 


(82 -5)/115 


1.31749(9) [17] 


2.07 bIAl,=i/4 


(35-3)/48 


1.31786(14) [30] 


bIAl,=i/2 


(28 - 2)/48 


1.31785(10)[29] 


bIAl,=i 


(19 - l)/48 


1.31784(10)[28] 


bIAlp=3^s=l/2 


(103 - 1)/172 


1.31766(21)[24] 


b±IAl^=i/4 


(36-2)/48 


1.31789(20)[28] 


b±IAl,=i/2 


(21 - l)/48 


1.31789(22)[28] 


bFCIAl,=i/4 


(36-4)/48 


1.31775(10)[28] 


bFCIAl,=i/2 


(35-4)/48 


1.31775(9) [28] 


bIA2,=i/8 


(99 - 7)/115 


1.31780(9) [27] 


bIA2,=i/4 


(93-4)/115 


1.31780(9) [27] 


bIA2,=i/2 


(87-4)/115 


1.31780(8) [28] 


bIA2s=i 


(60 - 2)/115 


1.31781(7) [27] 


bIA2„=i9 ,,=1/2 


(48-6)/70 


1.31777(10)[28] 


bIA2„=i8,s=i/2 


(53-4)/62 


1.31768(9) [28] 


bIA2p=3^s=i/2 


(277- 18)/345 


1.31772(14)[25] 


bIA2p=3,s=i 


(192 - 5)/345 


1.31773(14)[25] 


b±IA2,=i/2 


(46 - 3)/100 


1.31781(29)[23] 


bFCIA2,=i/2 


(91 - 2)/140 


1.31780(11)[29] 


bIA3,=i/2 


(56 - 4)/61 


1.31787(8) [31] 


2 10 hTAl 1 /n 


(r>q — '2) /AR 


1 31777('9'lf16l 


bIA2,=i/2 


(92 -2)/115 


1.31773(6) [15] 


bIA2p=3^s=i/2 


(295 - 17)/345 


1.31769(10)[14] 


b±IA2,=i/2 


(49 - 5)/100 


1.31774(20)[15] 


bFCIA2,=i/2 


(92 - 5)/140 


1.31772(15)[17] 


2.20 bIAl,=i/2 


(31 -3)/48 


1.31809(7) [30] 


bIA2,=i/2 


(94-6)/115 


1.31807(3) [27] 



From the intermediate results reported in Tables PCXIII| , PCXIV] , PCXV| and PCXVI] (which, 
we stress, are determined algorithmically once chosen the set of lA/c's), we obtain the esti- 
mates of 7 and u. 

From the analyses for the 0^ Hamiltonian at A = 2.07, we obtain 



7 = 1.31780(10) [28] + 0.003(A - 2.07), 
u = 0.67161(5)[12] + 0.002(A - 2.07). 



(B20) 
(B21) 



As before, the number between parentheses is basically the spread of the approximants at A = 
2.07 using the central value of f3c, while the number between brackets gives the systematic 
error due to the uncertainty on j3c- Eqs. ( |B2U| ) and ( [B21D show also the dependence of the 
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TABLE XXIV. Results for v obtained from the analysis of the 19th-ordcr HT series of ^ j fi 
for the (/>^ Hamiltonian. The number n of terms used in the analysis is indicated explicitly when 
it is smaller than the number of available terms {n = 19). p = when its value is not explicitly 
given. 



A 


approximants 




V 


2.00 


bIAl^=i/2 


36/37 


0.67140(2) [9] 




bIA2,=i/2 


(63-6)/70 


0.67141(4)[8] 


2.07 


bIAl,=i/2,i 


36/37 


0.67161(2)[13] 




bIAl„=ig ,,^1/2 


30/36 


0.67160(4) [13] 




bIAl„=17^s=l/2 


(30 - l)/33 


0.67163(11)[12] 




bIAlp=3_5=l/2 


(124-5)/134 


0.67162(5) [12] 




bIAlp=3,s=i 


(120-6)/134 


0.67162(5)[12] 




b±IAl,=i/2 


(33 - l)/36 


0.67161(2)[13] 




bFCIAl,=i/2 


(29 -3)/37 


0.67158(10)[12] 




bIA2,=i/2 


(66 - 5)/70 


0.67161(4) [12] 




bIA2s=i 


(50-3)/70 


0.67162(4) [12] 




bIA2„=i8,s=i/2 


(44-3)/62 


0.67162(5)[12] 




bIA2„=i7^s=i/2 


(38 - 2)/49 


0.67166(4) [11] 






(^ sn — (\\ /91 ^ 


n fi7i (\A-ii\\ \^ 91 

U.U I -LU^l U ) 1 J-^J 




bIA2p=3,^=l 


(145 - 7)/215 


0.67164(5) [12] 




b±IA2,=i/2 


(55 -3)/55 


0.67161(3) [13] 




bFCIA2,=i/2 


(60-4)/85 


0.67161(11)[14] 




bIA3,=i/2 


(17- l)/34 


0.67159(6) [14] 


2.10 


bIAl,=i/2 


36/37 


0.67160(2) [8] 




bIA2,=i/2 


(63-5)/70 


0.67161(4) [8] 


2.20 


bIAl,=i/2 


36/37 


0.67182(3) [14] 




bIA2,=i/2 


(60-4)/70 


0.67183(7) [14] 



results on the chosen vahic of A. The coefficient is estimated from the results for A = 2.2 
and A = 2.0, i.e., from the ratio 

g(A ^ 2.2) - g(A ^ 2.0) 

0.2 ' ^ ' 

where Q represents the quantity at hand. Using A* = 2.07(5), we obtain finally 

7 = 1.31780(10)[28]{15}, i/ = 0.67161(5)[12]{10}, (B23) 

where the error due to the uncertainty on A* is reported between braces. 

Since for A = 2.10 a more precise estimate of /?,. is available, it is interesting to perform 
the same analysis, using the HT series of the 0^ model at A = 2.10. We obtain 

7 = 1.31773(10) [15] + 0.003(A - 2.10), (B24) 
V = 0.67160(5)[8] + 0.002(A - 2.10), (B25) 
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TABLE XXV. 


Results for 7 obtained from the analysis of the 20th-order HT series of % 


for 


the dd-XY model. 


p = Q when its value is 


not explicitly given. 






D 


approximants 


ra 


7 




0.90 


bIAl,=i/2 


(35 - l)/48 


1.31685(29) 


[24] 




bIA2,=i/2 


(66 - 3)/115 


1.31693(28) 


[27] 


1.02 


bIAl,=i/4 


(46 - 2) /48 


1.31745(20) 


[21] 




bIAl,=i/2 


(41 - 3)/48 


1.31746(17) 


[22] 




bIAl,=i 


(24 - l)/48 


1.31748(15) 


[23] 




bIAlp=3 ,,=1/2 


(162 - 9)/172 


1.31733(35) 


[20] 




b±IAl,=i/2 


(35 - 2)/48 


1.31735(13) 


[22] 




bFCIAl,=i/2 


(31 - 4)/48 


1.31745(21) 


[22] 




bIA2,=i/4 


(103 - 3)/115 


1.31748(25) 


[23] 




bIA2,=i/2 


(68 - 1)/115 


1.31748(16) 


[22] 




bIA25=i 


(22 - 1)/115 


1.31754(20) 


[20] 




bIA2p=3^s=i/2 


(259 - 7) /345 


1.31730(26) 


[19] 




b±IA2,=i/2 


(74 - 3)/100 


1.31754(26) 


[22] 




bFCIA2,=i/2 


(69 - 4)/140 


1.31738(38) 


[18] 




bIA3,=i/2 


(41 - 1)/61 


1.31776(19) 


[24] 


1.03 


bIAl,=i/2 


(40 - 3)/48 


1.31747(15) 


[14] 




bIA2,=i/2 


(71 - 1)/115 


1.31751(13) 


[15] 




bIA2„ q „ 1 /n 


(263 - ll)/345 


1.31733(22) 


[13] 




b±IA2,=i/2 


(73-2)/100 


1.31756(24) 


f13l 




bFCIA2,=i/2 


(72 - 4)/140 


1.31744(27) 


[11] 




bIA3,=i/2 


(41 - 2)/61 


1.31776(16) 


[16] 


1.20 


bIAl,=i/2 


(43-l)/48 


1.31867(20) 


[28] 




bIA2,=i/2 


(99-4)/115 


1.31868(10) 


[25] 



which, using A* — 2.07(5), give 

7 = 1.31764(10)[15]{15}, i/ = 0.67154(5) [8]{10}, (B26) 

in perfect agreement with the results obtained at A = 2.07. The slight difference of the 
central values is essentially due to the independent estimates of (5c. 
Prom the analyses for the dd-XY model at = 1.02, we have 

7 = 1.31748(20) [22] + 0.006(L» - 1.02), (B27) 
u = 0.67145(10)[10] + 0.005(L' - 1.02), (B28) 

where the coefficient determining the dependence of the results on D is estimated by com- 
puting 

Q{D ^ 1.2) - Q(D ^ 0.9) 

0:3 ■ ^^^^^ 

Since D* — 1.02(3), we obtain the final estimates 



51 



TABLE XXVI. Results for v obtained from the analysis of the 19th-order HT series of 
for the dd-XY model, n = 19 and p = when not explicitly given. 



D 


approximants 


ra 


V 


0.90 


bIAl,,=i/2 


(33-2)/37 


0.67091(6) [12] 




bIA2,=i/2 


(62-3/70 


0.67092(10)[12l 


1.02 


bIAl,,=i/2 


(35 -3)/37 


0.67146(7) [101 

\ / L J 




bIAl,,=i 


(30 - l)/37 


0.67147(5) [101 

\ / L J 




n — io,s — -L/z 


(32 - l)/36 


0.67148(15)[10] 




blAl„_i7 o_i /9 


(31 - 2) 733 


0.67132(28 91 




bIAl„_o o_i /o 


(124-11)7134 


0.67143(12)[10l 




bIAlp=3^s=i 


(103 - 9)7134 


0.67145(10)fl0l 

V / L J 




b±IAl„_i/2 

^ 6 — J- / Z 


(34 - l)/36 


0.67143(5) [11] 




bFCIAl,=i/2 


(33 - 3)737 


0.67136(12)[10] 

V / L J 




bIA2,=i/2 


(64 - 1)770 


0.67144(5) [10] 




bIA2s=i 


(55-3)/70 


0.67145(4) [10] 




bIA2„_ic ,_i /2 


(54- 2)/62 


0.67145(11)[10] 




bIA2„=i7 s=x/2 


(48- 5)/49 


0.67137(3) [9] 




bIA2p=3 s=i/2 


(198 - 9)/215 


0.67141(8) [9] 




b±IA2,=i/2 


(53-9)/55 


0.67144(2) [10] 




bFCIA2,=i/2 


(78-8)/85 


0.67140(6) [11] 




bIA3,=i/2 


(34-4)/34 


0.67149(5) [10] 


1.03 


bIAl,=i/2 


(34-2)/37 


0.67149(7) [8] 




bIA2,=i/2 


(67-4)/70 


0.67147(5) [7] 


1.20 


bIAl,=i/2 


(30 - l)/37 


0.67231(12)[13] 




bIAl,=i/2 


(64-4)/70 


0.67236(7) [12] 




7 = 1.31748(20)[22]{18}, 


u = 0.67145(10)[10]{15}. 


(B30) 


Since for D = 


= 1.03 a more precise estimate of (3c is available, it is worthwhile to repeat the 


analysis usinj 


I the series at this value of D. We have 






7 = 1.31751(20)[15] + 0.006(L> - 1.03), 


(B31) 




u = 0.67148(10)[8 


] + 0.005(D- 1.03), 


(B32) 


which, for -D" 


' = 1.02(2), give 








7 = 1.31745(20)[15]{18}, 


z/ = 0.67143(10)[8]{15}, 


(B33) 



in good agreement with the results ( |B30| ). 

Consistent, although significantly less precise, results are obtained from IHT analyses 
that do not make use of the MC estimate of (5c- For example, by analyzing the HT series 
for the 0^ Hamiltonian at A = 2.07, we find /^^ = 0.509385(8), 7 = 1.3178(8){3}, v = 
0.6716(4){1}, where the error in parentheses is the spread of the approximants and the error 
between braces corresponds to the uncertainty on A*. Here, we determine (3c and 7 from the 
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analysis of using IA2's, FCIA2's, and IA3's, and v from the analysis of ^ using bIA2's 
biased with the estimate of /3c obtained in the HT analysis of x- 

From the results for 7 and z/, one can obtain 77 by the scaling relation 7 = (2 — r^)//. This 
gives ?7 = 0.0379(10), where the error is estimated by considering the errors on 7 and v as 
independent, which is of course not true. We can obtain an estimate of ri with a smaller, yet 
reliable, error using the so-called critical-point renormalization method (CPRM) (see Ref. 
10 and references therein). In the CPRM, given two series D{x) and E{x) that are singular 
at the same point Xq, D{x) = Ylii^i^^ ~ (^0 ~ and E{x) = ~ (^0 ~ ^j)"", one 

constructs a new series F{x) = ^.((ij/ej)x*. The function F{x) is singular at x = 1 and for 
X ^ 1 behaves as F{x) ~ (1 — x)'"^, where (p = 1 + 5 — e. Therefore, the difference 6 — e can 
be obtained by analyzing the expansion of F{x) by means of biased approximants with a 
singularity at Xc = 1. In order to check for possible systematic errors, we applied the CPRM 
to the series of f] and x (analyzing the corresponding 19th-order series) and to the series 
of and X (analyzing the corresponding 20th-order series). We used lA's biased at Xc = 1. 
In Table pCXVII| we present the results of several sets of lA's. For the 0^ model at A = 2.07 



we obtain 

r]i^ = 0.02550(20) + 0.0013(A - 2.07). (B34) 

Thus, taking into account that A* = 2.07(5), we find 

r]u = 0.02550(20){7}, (B35) 

where the first error is related to the spread of the lA's and the second one to the uncertainty 
on A*, evaluated as before. Analogously, for the dd-XY model we find 

r]u = 0.02550(40) + 0.004(L) - 1.02), (B36) 

and therefore, using D* = 1.02(3), 

r]u = 0.02550(40){12}, (B37) 

where again the first error is related to the spread of the lA's, while the second one is related 
to the uncertainty on D*. 



3. Amplitude ratios 

In the following we describe the analysis method we employed to evaluate zero- 
momentum renormalized couplings, such as and r2j. In the case of we analyzed the 
series p^^^g^ = YlZo^iP'- 

Consider an amplitude ratio A which, for t = Pc/ P — 1 — >■ 0, behaves as 

A(t) =A* + cit^ + C2t^^ + .... (B38) 

In order to determine A* from the HT series of A{t), we consider biased lAl's, whose behavior 
at Pc is given by (see, e.g., Ref. 0) 
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TABLE XXVII. Results for t] obtained using the CPRM: (a) applied to ^V/^ and x (19 orders); 
(b) applied to and x (20 orders). 





approximants 


ra 




A = 2.00 


(a) bIAl^=i/2 

(a) bIA2,=i/2 

(b) bIA2,=i/2 


33/37 
(47- l)/70 
(99 - 9)/115 


0.02547(7) 

0.0256(2) 

0.0251(3) 


A = 2.07 


(a) bIAl^=i/2 
(a) bIA2,=i/2 

(a) bIA3^=i/2 

(b) bIA2,=i/2 
(b) bIA3,=i/2 


37/37 

47/70 
(20 - l)/34 
(96 - 8)/115 
(51 - 2)/61 


0.02555(7) 

0.0257(2) 

0.0255(2) 

0.0252(3) 

0.0253(5) 


A = 2.20 


(a) bIAl^=i/2 

(a) bIA2,=i/2 

(b) bIA2,=i/2 


33/37 
(49 - 2)/70 
(95 - 11)/115 


0.02573(7) 

0.0259(3) 

0.0252(3) 


D = 0.90 


(a) bIA2,=i/2 

(b) bIA2,=i/2 


(45 - l)/70 
(84-3)/115 


0.0252(3) 
0.0248(9) 


n 1 no 

U — L.UZ 


(aj biAis=i/2 
(a) bIA2,=i/2 

(a) bIA3^=i/2 

(b) bIA2,=i/2 
(b) bIA3,=i/2 


too 1 \ /T? 

(37- l)/70 
(23 - 2)/34 
(93 - 3)/115 
(59 - 3)/61 


u.uzoD^y ) 
0.0257(3) 
0.0252(5) 
0.0252(8) 
0.0253(4) 


D = 1.20 


(a) bIA2,=i/2 

(b) bIA2,=i/2 


(33 - 2)/70 
(96 - 4)/115 


0.0263(3) 
0.0259(8) 



(B39) 



where /(/3) and g{P) are regular at /3c, except when ( is a non-negative integer. In particular, 



c 



PoiPc) 



gWc 



(B40) 



In the case we are considering, ( is positive and therefore, g{f3c) provides an estimate of A*. 
Moreover, for improved Hamiltonians we expect C = A2 ~ 2A and A ^ 0.5. In our analyses 
we consider blAl's and bj-IAl's (see Eqs. ( Bll ) and ( [B12|) ) and impose various constraints 
on the value of C by selecting blAl's with ( larger than a given nonnegative value. 

In Table pCXVllI| we report the results obtained for (74 using different sets of approximants. 
In this case the variation due to the uncertainty of (3c is negligible. Therefore, we report 
only the average of the results of the "good" lAl's and their standard deviation (divided 
by 



/Ta) calculated at Pc- In Table PCXVIII| we also report the value of C obtained from the 
selected lAl's. The comparison of the results for different values of A and D shows that the 
errors due to uncertainty on A* and D* are small and negligible. 



From the results of Table XXVIII we derive the estimates 



^4 = 21.15(6), 



^4 = 21.13(7), 



(B41) 
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TABLE XXVIII. Results for (74, obtained from the analysis of the 17th-order series of 







approximants 
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c 


X = 


2.00 




(41 - l)/43 


21.16(7) 


1.1(5) 






b±IAls=i/4^^>o 


(40 - 2)/44 


21.14(6) 


1.3(8) 


A = 


2.07 


t'IAl5_]^/xo,i/4,c>o 


(41 - 2)/43 


21.17(6) 


1.2(5) 






bIAl^ = ;^/4^^>0.5 


(38 - 2)/43 


21.16(6) 


1.2(5) 






bIAl^^;^/4 ]^ 3>^>Q 7 


23/43 


21.19(5) 


1.0(2) 






bIAlp=2,s=l/4,C>0 


(105 - 4)/118 


21.14(7) 


1.7(8) 






b±IAls=i/io,i/4,c>o 


(40-3)/44 


21.14(5) 


1.4(9) 






b±IAls=i/4,^>0.5 


(39 - 2)/44 


21.14(5) 


1.4(9) 






b±IAls=i/4^1.3>^>0.7 


(20- l)/44 


21.16(3) 


1.1(1) 






b±IAl„=2,s=l/4,C>0 


(80 - 3)/97 


21.13(7) 


2(2) 


A = 


2.10 


bIAls=i/4^^>o 


(41-2)/43 


21.17(6) 


1.2(5) 






bIAl5=l/4 ,j>0.5 


(38 - l)/43 


21.17(6) 


1.2(5) 






bIAls=i/4 1 3>^>0.7 


23/43 


21.19(5) 


1.0(2) 






bIAlp=2,s=l/4,C>0 


(105 - 4)/118 


21.14(7) 


1.7(8) 






b±IAls=i/4,^>o 


(39 - 2)/44 


21.16(6) 


1.4(6) 






b±IAlp=2,s=l/4,C>0 


(79 - 3)/97 


21.14(7) 


2(2) 


A = 


2.20 


bIAl^=i/4^^>o 


(40 - 2)/43 


21.19(5) 


1.3(6) 






b±IAls=i/4,^>o 


(39 - 2)/44 


21.17(5) 


1.5(9) 


D = 


0.90 




(32 - l)/43 


21.07(12) 


2(2) 






b±IAl5=i/4 ,j~^Q 


(35-5)/44 


21.07(9) 


1.1(4) 


D = 


1.02 


bIAl^=;^/io,c>o 


(31 - 2)/43 


21.16(10) 


1.5(1.4) 






bIAl5=;^/4^^>o 


(30-2)/43 


21.16(10) 


1.5(1.4) 






bIAl5=i/4^^>0.5 


(24- l)/43 


21.13(6) 


1.7(1.4) 






bIAl5=l/4 X.3>f>0.7 


9/43 


21.16(7) 


0.9(2) 






bIAlp=2,s=l/4,C>0 


(69 - 8)/118 


21.2(3) 


1.6(1.3) 






b±IAls=i/4,^>o 


(36-3)/44 


21.13(7) 


1.5(1.0) 






b±IAls=i/4,c>o.5 


(32 - 4)/44 


21.11(5) 


1.5(8) 






b±IAls=i/4,1.3>f>0.7 


(16- l)/44 


21.13(3) 


1.1(1) 






b±IAlp=2,s=l/4,C>0 


(66 - 7)/97 


21.13(12) 


2(2) 


D = 


1.03 


bIAls=i/4,^>o 


(32 - 3)/43 


21.17(12) 


1.5(1.5) 






b±IAl^=i/4^^>o 


(36-3)/44 


21.13(6) 


1.6(1.1) 


D = 


1.20 


bIAls=i/4,^>o 


(34-4)/43 


21.23(3) 


3(3) 






b±IAls=i/4,^>o 


(33 - 2)/44 


21.23(3) 


2(2) 



respectively for the 0^ Hamiltonian and the dd-XY model. We note that these results are 
slightly larger than the estimates reported in Ref. p6[. The difference is essentially due to the 



different analysis employed. There, the analysis was based on Fade (PA), Dlog-Pade (DPA) 
and lAl's, selecting those without singularities in a neighborhood of /3c and evaluating them 
at /?c- However, by analyzing the longer series that are now available for the Ising model 
|79|, we have realized that such procedure is not very accurate and that the analyses using 
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blAl's are more reliable when a sufficiently large number of terms is available. Moreover, 
when the series is sufficiently long, most (and eventually all) PA's, DPA's and lAl's become 
defective. Indeed, the functions we are considering do have singularities at /3c, although with 
a positive exponent. 

In the analysis of r2j, we also consider PA's and DPA's. We indeed expect that, when 
the series is not sufficiently long to be asymptotic, the approximants obtained by biasing 
the singularity at (3c may not provide a robust analysis. For comparison, we also use quasi- 
diagonal Pade approximants (PA's) and Dlog-Pade approximants (DPA's), evaluating them 
at Pc- 

We consider: 

[l/m] PA's with 

I + m > n — 2, 

Max[n/2-q,4]<l,m<n/2 + q, (B42) 

where I, m are the orders of the polynomials respectively in the numerator and denominator 
of the PA. As estimate from the PA's we take the average of the values at jSc of the nonde- 
fective approximants using all the available terms of the series and satisfying the condition 
with q = 3. The error we quote is the standard deviation (divided by ^/ra) of the re- 



sults from all the nondefective approximants listed above. We consider defective those PA's 
that have singularities in the rectangle defined in Eq. ( [B18| ) with: Xmin = 0.9, Xmax = 1.01, 
and i/max = 0.1 for and rs; Xmin = 0, Xmax = 1-1, and ?/max = 0.1 for riQ. 
[l/m] DPA's with 

I + m > n — 3, 

Max [(n - l)/2 - g, 4] < /, m < (n - l)/2 + q, (B43) 

where I, m are the orders of the polynomials respectively in the numerator and denominator 
of the PA of the series of its logarithmic derivative. We again fix g = 3. The estimate with 
the corresponding error is obtained as in the case of PA's. We consider defective those DPA's 
that have singularities in the rectangle ( |B18| ) with Xmin = 0, Xmax = 1-01, and ?/max = 0.1. 

For rg and rg the above PA's and DPA's give results substantially consistent with those 
of blAl's, see Tables pCXlX| and pCXX|. Therefore the systematic error due to the use of PA's 



and DPA's should be smaller than the final quoted errors. From Tables XXIX and XXX we 



derive the estimates reported in Table pClV| . For rio we obtain only very rough estimates 



using essentially PA's: rio = —13(7) from the 0^ Hamiltonian and rio = —11(14) from the 
dd-XF model. 



APPENDIX C: UNIVERSAL AMPLITUDE RATIOS 

We give here the definitions of the amplitude ratios that are used in the text. They are 
expressed in terms of the amplitudes derived from the singular behavior of the specific heat 

CH = A^\t\-\ (CI) 

the magnetic susceptibility in the high-temperature phase 
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TABLE XXIX. Results for vq from the analysis of the corresponding 17th-order HT series. 





approximants 




ra 


A = 2.00 


PA 


19/21 


1.945(20) 




DPA 


12/18 


1.959(27) 


A = 2.07 


PA 


19/21 


1.947(22) 




DPA 


12/18 


1.962(26) 




bIAls=i/io,c>o 


(21 - l)/43 


1.98(10) 




bIAl5=i/io,c>o.5 


(17- l)/43 


1.955(9) 




bIAlp=2,s=i/io,c>o 


(58-3)/118 


1.99(9) 




bIAlp=2,s=l/10,C>0.5 


(42 - 2)/118 


1.956(10) 


A = 2.10 


PA 


19/21 


1.948(23) 




DPA 


12/18 


1.963(26) 




bIAl^=i/io,c>o 


(21 - l)/43 


1.98(7) 




bIAls=i/io,c>o.5 


(18 - l)/43 


1.959(19) 


A = 2.20 


PA 


19/21 


1.951(30) 




DPA 


13/18 


1.967(26) 


D = 0.90 


PA 


21/21 


1.934(15) 




DPA 


16/18 


1.936(11) 


D = 1.02 


PA 


20/21 


1.945(13) 




DPA 


12/18 


1.952(13) 




bIAl^=i/io,c>o 


(28- l)/43 


1.948(11) 




bIAl^=i/io,c>o.5 


(27- l)/43 


1.947(8) 




bIAlp=2,s=i/io,c>o 


(69- 1)/118 


1.947(17) 




bIAlp=2,s=l/10,C>0.5 


(65- 1)/118 


1.948(16) 


D = 1.03 


PA 


20/21 


1.946(8) 




DPA 


12/18 


1.954(16) 


D = 1.20 


PA 


21/21 


1.963(6) 




DPA 


16/18 


1.963(12) 



X = (C2) 

the zero-momentum four-point connected correlation function in the high temperature phase 

X4 = ^c+r^-^^^ (C3) 

the second-moment correlation length in the high-temperature phase 

e = /+^-^ (C4) 

the spontaneous magnetization on the coexistence curve 

M^B\tf, (C5) 
and of the susceptibility along the critical isotherm 
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TABLE XXX. Results for from the analysis of the corresponding 16-th order HT series. 





approximants 


ra 


'^8 


A = 2.00 


PA 


17/18 


1.32(14) 




DPA 


14/21 


1.34(22) 


A = 2.07 


PA 


17/18 


1.36(14) 




DPA 


13/21 


1.37(15) 




1 X A -1 

biAls=i/io,^>0 


(10 — l)/33 


1.39(13) 




KT A 1 

DlAlp=2,s=l/10,C>0 


(^00 — ^j/yi 


1 "jfiM K\ 
i.oUyltD ) 


A = z.iU 




1 Q /I Q 

io/ io 


i.oo(iz ) 






1 Q /01 

io/2i 


i.oo(^iz J 


\ o on 


T3 A 


io/ io 


i.4i(iz j 






io/ Zl 


i.4Zl^i / J 


D = 0.90 


PA 


16/18 


1.34(18) 




DPA 


18/21 


1.28(10) 


D = 1.02 


PA 


17/18 


1.50(14) 




DPA 


15/21 


1.47(6) 




KT A 1 . „ J. „ 
DlAi5=i/io,(^>0 


(-\r: 1 WQQ 
t 10 — 1 J / 00 


1 .^01 IZ t 




bIAlj,=2,s=i/io,c>o 


(45 - 5)/91 


1.50(33) 


D = 1.03 


PA 


17/18 


1.52(13) 




DPA 


15/21 


1.54 8) 


D = 1.20 


PA 


16/18 


1.64(9) 




DPA 


15/21 


1.65(13) 




XL = C^\H\-^/<^\ 




(C6) 


We consider the following universal amplitude ratios: 








aA+C+ 




(C7) 




^ - {C+f ~ ' °' ' 




(C8) 








(C9) 



R+ = (yl+) = (^^) . (CIO) 



58 



REFERENCES 



[1] J. A. Lipa, D. R. Swanson, J. A. Nissen, T. C. P. Chui, and U. E. Israelsson, Phys. Rev. 

Lett. 76, 944 (1996). 
[2] R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 
[3] M. Hasenbusch and T. Tor5k, J. Phys. A 32, 6361 (1999). 

[4] J. A. Lipa, D. R. Swanson, J. A. Nissen, Z. K. Geng, P. R. Wilhamson, D. A. Strieker, 

T. C. P. Chui, U. E. Israelsson, and M. Larson, Phys. Rev. Lett. 84, 4894 (2000). 
[5] J. A. Lipa, private communication. 

[6] M. Campostrini, A. Pehssetto, P. Rossi, and E. Vicari, Phys. Rev. B 61 5905 (2000). 
[7] J. A. Nissen, D. R. Swanson, Z. K. Geng, V. Dohm, U. E. Israelsson, M. J. DiPirro, 

and J. A. Lipa, Low Temp. Phys. 24, 86 (1998). 
[8] K. E. Newman and E. K. Riedel, Phys. Rev. B 30, 6615 (1984). 

[9] D. L. Hunter and G. A. Baker, Jr., Phys. Rev. B 7, 3346 (1973); B 7, 3377 (1973); B 
19, 3808 (1979); M. E. Fisher and H. Au-Yang, J. Phys. A 12, 1677 (1979); Erratum A 
13, 1517 (1980); A. J. Guttmann and G. S. Joyce, J. Phys. A 5, L81 (1972); J. J. Rehr, 
A. J. Guttmann, and G. S Joyce, J. Phys. A 13, 1587 (1980). 

[10] A. J. Guttmann, in Phase Transitions and Critical Phenomena, Vol. 13, edited by 
C. Domb and J. Lebowitz (Academic, New York, 1989). 

[11] R. Z. Roskies, Phys. Rev. B 24, 5305 (1981). 

[12] J. Adler, M. Moshe, and V. Privman, Phys. Rev. B 26, 1411 (1982); Phys. Rev. B 26, 
3958 (1982). 

[13] P. Butera and M. Comi, Phys. Rev. B 56, 8212 (1997). 

[14] A. Pelissetto and E. Vicari, Nucl. Phys. B 519, 626 (1998); Nucl. Phys. B (Proc. Suppl.) 
73, 775 (1999). 

[15] P. Butera and M. Comi, Phys. Rev. B 58, 11552 (1998). 

[16] M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. E 60, 3526 (1999). 

[17] J.-H. Chen, M. E. Fisher, and B. G. Nickel, Phys. Rev. Lett. 48, 630 (1982). 

[18] M. J. George and J. J. Rehr, Phys. Rev. Lett. 53, 2063 (1984). 

[19] M. E. Fisher and J.-H. Chen, J. Physique 46, 1645 (1985). 

[20] B. G. Nickel and J. J. Rehr, J. Stat. Phys. 61, 1 (1990). 

[21] M. Hasenbusch, K. Finn, and S. Vinti, Phys. Rev. B 59, 11471 (1999). 

[22] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, and A. Munoz Sudupe, Phys. 

Lett. B 441, 330 (1998). 
[23] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. Munoz Sudupe, G. Parisi, and 

J. J. Ruiz-Lorenzo, J. Phys. A 32, 1 (1999). 
[24] M. Hasenbusch, J. Phys. A 32, 4851 (1999). 
[25] M. Hasenbusch, in preparation. 

[26] M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 62, 5843 (2000). 

[27] P. H. Damgaard and M. Hasenbusch, Phys. Lett. B 331, 440 (1994). 

[28] F. Jasch and H. Kleinert, "Fast convergent resummation algorithm and critical expo- 
nents of 0^ theory in three dimensions," e-print |cond-mat / 990624B| v2. 

[29] G. A. Baker, Jr., B. G. Nickel, M. S. Green, and D. I. Meiron, Phys. Rev. Lett. 36, 
1351 (1977); G. A. Baker, Jr., B. G. Nickel, and D. I. Meiron, Phys. Rev. B 17, 1365 
(1978). 



59 



D. B. Murray and B. G. Nickel, "Revised estimates for critical exponents for the con- 
tinuum n- vector model in 3 dimensions," unpublished Guelph University report (1991). 
K. G. Chetyrkin, S. G. Gorishny, S. A. Larin, and F. V. Tkachov, Phys. Lett. B 132, 
351 (1983). 

H. Kleinert, J. Neu, V. Schulte-Frohlinde, K. G. Chetyrkin, and S. A. Larin, Phys. Lett. 
B 272, 39 (1991); Erratum B 319, 545 (1993). 
P. Butera and M. Comi, Phys. Rev. B 60, 6749 (1999). 
M. Krech and D. P. Landau, Phys. Rev. B 60, 3375 (1999). 

H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, and A. Munoz Sudupe, Phys. 
Lett. B 387, 125 (1996). 

N. Schuhka and E. Manousakis, Phys. Rev. B 52, 7528 (1995). 
M. Hasenbusch and A. P. Gottlob, Physica A 201, 593 (1993). 
W. Janke, Phys. Lett. A 148, 306 (1990). 

L. S. Goldner, N. Mulders, and G. Ahlers, J. Low Temp. Phys. 93, 131 (1993). 

D. R. Swanson, T. C. P. Chui, and J. A. Lipa, Phys. Rev. B 46, 9043 (1992); D. Marek, 
J. A. Lipa, and D. Philips, Phys. Rev. B 38, 4465 (1988). 
A. Singsaas and G. Ahlers, Phys. Rev. B 30, 5103 (1984). 
J. A. Lipa and T. C. P. Chui, Phys. Rev. Lett. 51, 2291 (1983). 
M. Hasenbusch, Physica A 197, 423 (1993). 
A. P. Gottlob and M. Hasenbusch, J. Stat. Phys. 77, 919 (1994). 

M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Europhys. Lett. 38, 577 (1997); 
Phys. Rev. E 57, 184 (1998). 

V. Privman, P. C. Hohenberg, and A. Aharony, in Phase Transitions and Critical Phe- 
nomena, Vol. 14, edited by C. Domb and J. L. Lebowitz (Academic, New York, 1991). 
R. Guida and J. Zinn- Justin, Nucl. Phys. B 489, 626 (1997). 

J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, third edition (Claren- 
don Press, Oxford, 1996). 

E. Brezin and D. J. Wallace, Phys. Rev. B 7, 1967 (1973). 
E. Brezin and J. Zinn- Justin, Phys. Rev. B 14, 3110 (1976). 

I. D. Lawrie, J. Phys. A 14, 2489 (1981). 
D. J. Wallace and R. P. K. Zia, Phys. Rev. B 12, 5340 (1975). 
L. Schafer and H. Horner, Z. Phys. B 29, 251 (1978). 
A. Pehssetto and E. Vicari, Nucl. Phys. B 540, 639 (1999). 
A. Pelissetto and E. Vicari, Nucl. Phys. B 575, 579 (2000). 
T. Reisz, Phys. Lett. B 360, 77 (1995). 

A. 1. Sokolov, E. V. Orlov, V. A. Ul'kov, and S. S. Kastanov, Phys. Rev. E 60, 1344 
(1999). 

A. Pehssetto and E. Vicari, Nucl. Phys. B 522, 605 (1998). 
P. Schofield, Phys. Rev. Lett. 22, 606 (1969). 

P. Schofield, J. D. Lister, and J. T. Ho, Phys. Rev. Lett. 23, 1098 (1969). 

B. D. Josephson, J. Phys. C: Sohd State Phys. 2, 1113 (1969). 

J. Engels, S. Holtmann, T. Mendes, and T. Schulze, "Equation of state and Goldstone- 
mode effects of the three-dimensional 0(2) model," e-print |hep-lat / 0006023] , Phys. Lett. 
B, in press. 

[63] C. Hohenberg, A. Aharony, B. I. Halperin, and E. D. Siggia, Phys. Rev. B 13, 2986 



60 



(1976). 

S. A. Larin, M. Monnigman, M. Strosser, and V. Dohm, Phys. Rev. B 58, 3394 (1998). 
C. Bervillier, Phys. Rev. B 34, 8141 (1986). 

T. Takada and T. Watanabe, J. Low Temp. Phys. 49, 435 (1982). 

M. E. Fisher and R. J. Burford, Phys. Rev. 156, 583 (1967); H. B. Tarko and 

M. E. Fisher, Phys. Rev. Lett. 31, 926 (1973); H. B. Tarko and M. E. Fisher, Phys. 

Rev. B 11, 1217 (1975); M. E. Fisher and A. Aharony, Phys. Rev. B 7, 2818 (1974). 

R. A. Ferrel and D. J. Scalapino, Phys. Rev. Lett. 34, 200 (1975). 

A. J. Bray, Phys. Rev. B 14, 1248 (1976). 

M. E. Fisher and J. S. Langer, Phys. Rev. Lett. 20, 665 (1968). 
U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

S. Caracciolo, R. G. Edwards, A. Pehssetto, and A. D. Sokal, NucL Phys. B (Proc. 

SuppL) 26, 595 (1992); NucL Phys. B 403, 475 (1993). 

R. C. Brower and P. Tamayo, Phys. Rev. Lett. 62, 1087 (1989). 

M. Hasenbusch, Habihtationsschrift, December 1999, Humboldt- Univer sit at zu Berlin. 
M. Caselle, F. Gliozzi and S. Necco, "Thermal operators and cluster topology in g-state 
Potts Models," e-print |cond-mat/ 0006204 



F. Gliozzi and A. Sokal, private communication. 

H. G. Ballesteros and V. Martin-Mayor, Phys. Rev. E 58, 6787 (1998). 
M. E. Fisher, Philos. Mag. 7, 1731 (1962). 

M. Campostrini, "Linked- Cluster Expansion of the Ising Model" , e-print 
|cond-mat/000513q . 



61 



